Quantized Response and Topology of Insulators with Inversion Symmetry 
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We study three dimensional insulators with inversion symmetry, in which other point group 
symmetries, such as time reversal, are generically absent. Their band topology is found to be 
classified by the parities of occupied states at time reversal invariant momenta (TRIM parities), 
and by three Chern numbers. The TRIM parities of any insulator must satisfy a constraint: their 
product must be +1. The TRIM parities also constrain the Chern numbers modulo two. When the 
Chern numbers vanish, a magneto-electric response parameterized by is defined and is quantized 
to 6 = 0, 7r. Its value is entirely determined by the TRIM parities. These results may be useful in 
the search for magnetic topological insulators with large 6. A classification of inversion symmetric 
insulators is also given for general dimensions. An alternate geometrical derivation of our results is 
obtained by using the entanglement spectrum of the ground state wave-function. 



I. INTRODUCTION 

Recently, the discovery of a class of insulator with 
nontrivial band topology protected by time reversal 
symmetry^ has led to a range of discoveries on the quali- 
tative properties of solids. Soon after, it was realized that 
this is part of a broader classification of gapped topolog- 
ical phases with free fermions. If one includes particle 
hole symmetry, natural for superconducting Hamiltoni- 
ans, one ends up with ten symmetry classes^ that in- 
clude both insulators and gapped superconductors. The 
symmetry transformations, time reversal and particle 
hole symmetry, can be consistently generalized to in- 
clude disorder and inhomogencitics. Topological phases 
in these classes are always characterized by protected sur- 
face states. However, one can also discuss band topology 
in the presence of the point group symmetry of a crystal. 
That is, which classes of Hamiltonians can be smoothly 
connected to one another while preserving the insulating 
gap and the symmetry? Although these distinctions re- 
quire considering a perfect crystal, and inevitably fade 
in the presence of disorder, one often deals with crys- 
tals that are sufficiently clean for such distinctions to be 
useful. 

The general classification is a daunting task, given the 
vast number of symmetry groups of crystals. For ex- 
ample, there are 230 space groups; if one adds transfor- 
mation under time reversal symmetry to describe mag- 
netic insulators, one has 1,651 different groups. We focus 
here on a simple case, inversion (r — > — r), a symmetry 
that is commonly realized in magnetic insulators. For ex- 
ample, all Bravais lattices are inversion symmetric. An 
additional virtue is that inversion symmetry forces the 
magnetoclectric response to be isotropic and quantized, 
allowing us to calculate it by qualitative methods. Here, 
we will classify the topologically distinct inversion sym- 
metric insulators, and determine their protected proper- 
ties. (Phases of antiferromagnets can also be classified 
using a similar approach^.) 

The classes of three dimensional inversion symmetric 
insulators are found to be parameterized by three Chern 
numbers and a set of inversion parities. The Chern num- 



bers, which are already present in the absence of inversion 
symmetry, determine the quantized Hall response. While 
this is an integer in two dimensions, it is given by a re- 
ciprocal lattice vector Gh in a three dimensional crys- 
tal. Inversion parities are the only additional parameters 
that appear from adding the extra symmetry. They are 
defined at those momenta that are left invariant under 
k — > —k, up to a reciprocal lattice vector. This defines 
eight points in the 3D Brillouin Zone - the Time Rever- 
sal Invariant Momenta (TRIM). States at these momenta 
can be classified by their parities, eigenvalues under in- 
version, which are ±1. The set of inversion parities of the 
filled levels at the TRIMs and Gh completely specifics 
the topological class of the insulator. 

Usually, distinct topological classes can be distin- 
guished by their surface state content. A different ap- 
proach is to identify a quantized response, such as Hall 
conductivity, that can distinguish topological classes. For 
inversion symmetric insulators, in general, no protected 
surface states arise. This is due to the fact that although 
inversion is present in the bulk, it is absent at the sur- 
face of the system, which separates an inside from an 
outside. However, quantized physical response functions 
exist, and are protected by inversion symmetry. These 
include the magnetoelectric susceptibility 9, which is de- 
fined, for example, as the polarization (P) induced by a 
magnetic field (B), in a parallel direction: P = 
The coefficient 9 requires both bulk and surface to be 
insulating and is defined modulo 2ir. Under time rever- 
sal, 9 — > —9. Apart from the trivial solution 9 = 0, the 
ambiguity in the definition of 9 allows also for 9 = tt. 
The latter was shown to occur in time reversal symmet- 
ric topological insulators, when magnetic perturbations 
restricted to the surface open a surface gap£~— . 

Under inversion symmetry, 9 transforms in the same 
way as under time reversal. Once again we expect 
9 = 0, tt. The latter case is particularly interesting since 
the magnetoelectric susceptibility is large and no extra 
effort will usually be required to gap out surface states, 
in contrast to topological insulators. Typically, one does 
not expect surface modes in the inversion symmetric case. 
The response 9 is determined just by the parities of the 
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states at TRIMs, according to a formula derived below 
(which is analogous to formulae for the time reversal- 
symmetric cases&&2). Thus it will be possible to de- 
termine the magnetoelectric susceptibility without doing 
any band integrals. This should be useful in the search 
for large 9 insulators in magnetically ordered materials. 
Moreover, in magnetic insulators a large spin orbit cou- 
pling is not a priori required to obtain the phase with 
9 = 7r, which opens up a wider selection of candidates. 

A second set of results are in the form of constraints 
on the allowed values of TRIM parities and Chern num- 
bers Gh- We first show that in any insulator, a parity 
constraint is present. The product of the TRIM pari- 
ties of filled states is always +1. While this is a trivial 
condition when time reversal is also present, and energy 
levels come in Kramers pairs, it is rather non-trivial for 
magnetic insulators. For example, the parity assignment 
in Fig. [T]a, rules out a band gap. A corollary of this 
result is that it is impossible to make a direct continu- 
ous transition between a 9 = and 9 = tt insulator, if 
time reversal is broken. One must either evolve through 
a non-insulating state, or a first order transition. A con- 
tinuous transition only exists if both time reversal and 
inversion are present^. Another constraint fixes whether 
the three components of G^, are even or odd in relation 
to the TRIM parities. The intrinsic polarization P of the 
insulator is also determined by the parities. 

Finally we consider the entanglement. The entangle- 
ment spectrum of topological insulators contains special 
modesiir— . Section IIVI relates the number of gaplcss 
modes in the entanglement spectrum to the parity in- 
variants. This result is used to rederive some of the elec- 
tromagnetic properties in a simple fashion. 

The present article addresses some questions left open 
in earlier work in which we participated. Ref. [TH dis- 
cusses the entanglement spectrum of inversion symmet- 
ric insulators, without presenting the exact relation to 
TRIM parities. The proofs given here also complete 
the reasoning presented ir>i^ which studied the electronic 
structure of a specific material using the expressions for 
9 in terms of TRIM parities and the parity constraint. 

In the next section, we give a brief survey of all of our 
results and a guide to the sections where their proofs can 
be found. 



II. SUMMARY OF RESULTS 

The phases associated with a given symmetry can be 
classified using some topology. Each state can be classi- 
fied by integer-valued parameters (or by integers relative 
some modulus). These cannot change continuously, so 
they determine phases. 

Let us first define some conventions about the crystal 
lattice. We will for simplicity assume that the lattice is 
cubic (although there is no symmetry beyond inversion) 
and has a lattice spacing equal to one unit. All quanti- 
ties will be written with respect to a coordinate system 



xyz that is aligned with the axes of the crystal. The 
results may be translated to systems with a general Bra- 
vais lattice, by interpreting the expressions in the right 
coordinate system. Let be the primitive vectors of the 
lattice and let g* be the reciprocal vectors, g l -Rj = 2^5* . 
If a vector is in real space, the coordinates v x ,v y ,v z re- 
fer to v = v x Ri + v v R2 + v z 'R,3. Vectors u in reciprocal 
space should be expanded in terms of reciprocal vectors, 
u = j^iuxg 1 + u y g 2 + u z g 3 ). Electrical polarization is a 
real space vector; 3D Hall conductivity and momentum 
are in reciprocal space (and we also use upper and lower 
indices on the coordinates as a reminder of what basis to 
use). 

To classify bulk insulators, it is useful to look at inver- 
sion parities, as in the study of spectra of small molecules. 
The main difference is that in solids, the occupied states 
can be labelled by momentum. Let these states be given 
by ^ik(r) = Wik(r)e r . The parity classification is only 
useful at the "TRIMs," the momenta given by 

ni n 2 n 3 
K = ygi + y §2 + ygs (1) 

where ra!, 712,713 are integers. Such a momentum maps 
to itself under inversion symmetry modulo the reciprocal 
lattice, — k = k. Hence the wave functions at k must be 
invariant, and their parities can be defined: 

T-^an (r) = T) a («) V>a« (t) (2) 

Appendix [B] explains how to find these parities using a 
tight-binding model. 

We now introduce a key quantity u (k) at every TRIM 
k. This is defined as the number of states with odd par- 
ities at that TRIM. Note, these cannot change without 
a phase transition (at least in a non- interacting system) . 
Besides these 8 integers, the quantum Hall conductance 
gives three more invariant integers, since it is quantized: 
Gh = -§xh£*H where ^f- has integer components (ac- 
cording to the conventions defined above). 

These 11 integers, together with the total number of 
occupied bands n, are the only parameters necessary 
to determine a phase-any two band structures with the 
same integers can be tuned into one another without a 
phase transition. This scheme is derived in Sec. MI Al 
Appendix [A] gives an alternative method that is easier 
to generalize. These integers cannot be chosen indepen- 
dently of one another; there are some relationships be- 
tween their parities. These relationships can help predict 
material properties, since the parities are easy to deter- 
mine from the band structure. 

Total Parity Constraint:. All the constraints can be 
written in terms of the net parities: 

^ = (-ir (K) = ( 3 ) 

a 

For any insulator, one can show: 
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FIG. 1: Determining properties of systems using parities. The 
boxes represent an eighth of the Brillouin Zone; the TRIMs 
are at the corners. The signs represent the parities of the 
occupied states at the TRIMs. In (a) the parity constraint of 
even number of odd parity states is violated, hence it cannot 
be an insulator. In (b) the parities require a nonvanishing Hall 
conductance, with odd Chern number in the k x k y planes, (c) 
Quantized magnetoelectric response 9 — tv determined from 
number of odd parity states being 2 (mod 4). (d) A parity 
configuration corresponding to a frozen polarization. 

That is, the total number of filled odd parity states must 
be even. This is shown in Sec. MI Bl 

If a system has the parities in Fig. QJi, it must be 
metallic, because the parities do not satisfy Eq. (QJ. The 
gap must close at some momentum k in the Brillouin 
zone. The simplest case is if the system is a semimetal, 
with a single pair of cone-points at the Fermi energy. 

Quantum Hall Effect. Furthermore, the parities deter- 
mine the quantum Hall integers modulo 2. The parity of 
%pi is constrained; e.g. the z-component satisfies 

eiGnz = ]J Vk (5) 

k-R z =0 

That is, whether the Hall conductivity along the in- 
direction is an even or odd multiple of 2-k can be deter- 
mined by multiplying the 17's around either of the squares 
parallel to the xy-plane. This result is derived in Sec. 

ITncl 

If a system has the parities shown in Fig. [T}}, the Hall 
conductivity cannot vanish. The component along the z- 
direction, Ghz , must be an odd multiple of #- (per layer 
of the crystal). 

Eqs. (HJ),© are the only relationships between the 



invariants-if the relationships arc satisfied, the invariants 
can be realized in principle. 

Magnetoelectric Effect. The most interesting electro- 
magnetic response that we discuss is the magnetoelectric 
polarizability Oj. An applied magnetic field induces a 
polarization, P l = ctjB^ . In the absence of the quantum 
Hall effect, a* is well-defined. (Otherwise the polariza- 
tion can be neutralized by a flow of charge in the surface 
states associated with the Hall effect.) The polarizability 
a l j is odd under inversion: Under inversion symmetry, P 
changes sign, while B does not. 

If the crystal is inversion symmetric, it seems that a 
must vanish. However, a is ambiguous like the polariza- 
tion. An isotropic portion (j-8ij x integer) is indetermi- 
nate. (The magnetoelectric effect affects only properties 
of the surface of a crystal. An integer magnetoelectric 
effect can be mimicked by a quantum Hall coating on 
the surface.) Thus a* can be isotropic, ^j^9Sp if is a 
multiple of tt. To determine 9, the rj K parameters are not 
sufficient, and we must return to the n 's: 

- = 5$>°(") ( m °d 2). (6) 

According to Eq. ((4]), this is always an integer. This 
expression is proved in Sec. MI Dl 

Fig. [Tfc shows the parities of a model with 2 filled 
bands, which has a nontrivial 9. The total number of 
odd states at all the TRIMs is twice an odd number, 
hence 9 is 7T. One interpretation of 9 — tt is that the 
insulator has a half-integer quantum-Hall effect on the 
surface. Two filled bands are necessary for a nontrivial 
9. The requirement G# = implies, for one filled band, 
that the number of states is a multiple of 4. (See Eq. ([7] 
in the next section.) 

Frozen Polarization. Finally, let us discuss what in- 
formation is contained in the "net parities" r/ K . Eq. ([5]) 
shows that they determine the Chern numbers modulo 2. 
This accounts for three of the eight parities. Note that 
any pattern of parities satisfying the constraint Eq. (|4]) 
can be factored into 7 basis patterns: 

1](k) = i(-l)2?' G,, '' ! " s - + ^ i ' !, ' i *^^ H "' i *' I| ''(-l)»^ e '' ! 

(7) 

where the components of P e arc half-integers and the 
components of Gh arc integers times 2ir. (The factors of 
7T are added to make the exponents into integers.) This 
equation agrees with Eq. so Gh is the Hall con- 

ductivity modulo 2. When Gh = 0,^ varies as a plane 
wave on the vertices of the cube. The components of P e 
determine the intrinsic polarization. 

A crystal may have an intrinsic polarization when the 
Hall conductivity is zeroi^; in this case, r/^ looks like 
a plane wave on the vertices of the cube (see Fig. [TJi) 
and P e is the wave-number. The spontaneous electrical 
polarization is similar to 9; it is defined onlyi^ modulo 
a lattice vector times e. Inversion symmetry constrains 
the components to be integers or half-integers times e. 



Hence the polarization is determined by three bits and 
Eq. ([7J gives all the information about the polarization 
that can be obtained from bulk properties. To get the 
polarization correct, one needs to include the offset be- 
tween the electrons and the compensating charges in the 
nuclei: 



er Nl 



(8) 



where r^i is the position vector of the i th nucleus, with 
charge —Zie. This result is derived in Sec. IIII CI 

Consider the polarization of the crystal with the band 
structure illustrated in Fig. QJl; it is |Ri if the nuclei are 
all on the sites of the Bravais lattice. This intrinsic dipolc 
moment may correspond to actual ferroelectricity; the 
crystal would have a surface charge and a large electric 
field. Alternatively, the translational symmetry of the 
surface may be spontaneously broken or the surface may 
be metallic (see Ref. [l6| which is summarized in appendix 

IB- 

Parity constraints in general dimensions. The results 
in higher dimensions have a surprising feature: as the 
number of dimensions increases the sum of the n 's must 
be divisible by larger and larger powers of 2. 

Specifically, in 2s-dimensions, the sum of the n D 's is 
a multiple of 2 S_1 . This multiple is related to the 2s- 
dimensional Chern number G2 S (defined as a multiple of 
2tt): 



— y 

TRIM k 



u (k) 



s th Chern number 
2^ 



(mod 2); (9) 



the quantum Hall conductance is the 2-dimensional spe- 
cial case. 

In 2s + 1-dimensions, the sum of the n D 's is a multiple 
of 2 s and is related to the Chcrn-Simons integral 
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^^2n (K) (mod 2). 



(10) 



where the polarization and magnctoelcctric effect are the 
one- and three-dimensional versions. 

Note that insulators with inversion symmetry are quite 
different from ones: There is an insulator in 2s dimen- 
sions with a Chern number Gis equal to 1 which has 
just s filled bands2£ This insulator is not inversion sym- 
metric, though. The simplest inversion symmetric insu- 
lator with the identical Chern number has a minimum 
of 2 s bands, exponentially more bands than are necesary 
without symmetry. 

The Entanglement Spectrum. An insulator with in- 
version symmetry has a particle-hole symmetry I e in its 
entanglement spectrum e a (k) when it is cut on a plane 
through a center of inversion. This makes it very easy to 
determine qualitative properties of the Fermi arcs of the 
entanglement spectrum-it is possible to count (without 
topological arguments) the number of zero modes in the 




b) 




FIG. 2: Entanglement spectrum of a hopping Hamiltonian. 
a) The parities at the TRIMs. b) The entanglement modes 
on a cut parallel to the xy-plane. Note that there are two 
zero-modes at the TRIM (0, 0) and none at the other TRIMs, 
as expected from the parities. 



entanglement spectrum at the TRIMs k± along the sur- 
face. Let AN e (K±) = tr e= oX e ; that is AN e (nj_) is the 
number of even modes minus the number of odd modes 
with zero entanglement-energy at k±. 

This can be expressed in terms of parities of the bulk 
states through h (k). Define these parities relative to an 
inversion center on the plane of the cut. The quantity 
that appears is AJV(k) = tiE<oI, or n — 2n (n): 



AN e (K A 



-(AN( Ki ) + AN(k 2 )), 



(11) 



where K\ and k,^ are the two TRIMs that project to 
K±. In words: the difference between the number of even 
and odd states on the entanglement "Fermi surface" at 
a TRIM is half the difference between the even and odd 
states in the bulk, at the corresponding TRIMs. 

To illustrate an actual entanglement spectrum, we con- 
structed a Hamiltonian with a cubic unit cell whose in- 
version parities suggest that 8 = n and Ghi = (mod 2). 
The parities and the spectrum are shown in Fig. [5J (The 
Hamiltonian is described in Appendix |Bj) The entan- 
glement spectrum was calculated for a cut along the xy- 
plane. As expected, there is a Dirac point at (0,0). In 
this case, there are no chiral modes. 

From the relation between the entanglement spectrum 
and the parities, one can give alternative derivations of 
some of the results described above. This formula also 
leads to a simple determination of Fu and Kane's formula 
for the indices of topological insulators. These indices de- 
scribe the number (modulo 2) of physical surface states 
on a line connecting two TRIMs, and these states are 
easy to count for the entanglement spectrum. The en- 
tanglement spectrum can be continuously deformed into 
the physical spectrum so it has the same index. 
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III. ELECTROMAGNETIC EVALUATIONS 
AND EQUIVALENCE OF INSULATORS 

A. Classifying Inversion Symmetric Insulators 

This section presents the classification of noninteract- 
ing insulators with inversion symmetry. Hamiltonians 
with inversion symmetry can be related to Hamiltonians 
without any symmetry, which are already classified by 
Chern numbers. The only new parameters occurring for 
inversion symmetry come from the classification of " zero- 
dimensional insulators", or finite molecules with inver- 
sion symmetry. 

Consider the Hamiltonian H(k) for the wavefunctions 
VVik- This Hamiltonian can be taken to be an TV x N 
matrix by using a tight-binding model with N bands (n 
of which are filled). Since this Hamiltonian is inversion 
symmetric, 

I H(k)I = ff(-k) (12) 

where Iq is the matrix describing how the orbitals within 
the unit cell transform under inversion^ 7 -. 

Assume that H varies continuously without passing 
through a phase transition. This means that no states 
ever cross the Fermi energy, \i = 0, say. Let us determine 
when two Hamiltonians arc in the same phase-i.e., can 
be connected in this way. We wish to find criteria on the 
matrix fields H(k) that can be used to determine which 
phase they are in. 

The only special points in the Brillouin zones are the 
TRIMs, k, which are invariant under inversion symme- 
try. Each of these points can be interpreted as a zero- 
dimensional system, with a Hamiltonian H(k) that is 
invariant under Iq, since Eq. (|12[) implies IqH(k)Iq = 
H(k). Let n (K) be the number of eigenvalues at nega- 
tive energy which are odd under Iq. The states at this 
TRIM can mix together, but even states can mix only 
with even states and odd ones can mix only with odd 
ones, so the value of ti {k) cannot change. 

The Chern numbers of the Hamiltonian are topological 
winding numbers which also turn out to describe the Hall 
conductivit y 17 ' 18 . Because they are integers, they are 
also invariant. 

We will now spend the rest of this section showing that 
these give a complete classification of Hamiltonians with 
inversion symmetry. That is, if H(k) and H'(k) are two 
Hamiltonians with the same number n of occupied states, 
and the same total number of states N, and total number 
of odd states N a (both filled and empty bands) and such 
that 

ti (k) = n' {n) (for all TRIMs k) 



then (at least if N—n, n>2) there is a family of Hamilto- 
nians that connects H (k) to -ff'(k) without a phase tran- 
sition and while satisfying Eq. (fT2"|) . (We may assume 
that Iq remains constant since N, N Q do not change.) 



We do not usually consider the integers N and iV to 
be important invariants-their values can be changed by 
adding even or odd orbitals with a very high energy.In 
continuous space, there are infinitely many available or- 
bitals. 

The assumption N—n,n > 2 is included because, when 
there are too few bands, there are some Hamiltonians 
that cannot be deformed into one another just because 
there are not enough degrees of freedomi 1928 Our classi- 
fication theorem does not capture these distinctions, but 
the distinctions are not related to any generic properties. 
If one adds sufficiently many trivial occupied and unoc- 
cupied bands to an insulator, any two insulators with the 
same invariants can be deformed into one another. 

Proof for (|13p : Here we present the essential ideas 
in the classification of insulators. App. [5] gives a more 
systematic way of deriving the classification, including 
higher dimensions. 

This result can be derived by relating a Hamiltonian 
in d dimensions to one in a smaller number of dimen- 
sions. Let us take d to be arbitrary at first, so that 
we can describe the general reasoning. Let Hd be the 
space of general Hamiltonians in d dimensions, while Id 
is the subspace of Hamiltonians that also have inver- 
sion symmetry A generic Hamiltonian in %d can be 
regarded as a closed loop in Hd-i- For each value of kd 
(the d th component of k) , consider the d — 1-dimensional 
Hamiltonian Hk d defined by fixing one component of k, 
Hk d (ki, k 2 , ■ ■ ■ , kd-i) = H(ki, k 2 ,.--, k d ). This describes 
a closed loop because the Brillouin zone is periodic. 

A Hamiltonian in Id is an arc in Hd-i with end-points 
in Id-i (Fig. [3]) shows T-Ld schematically, with several 
possible arcs representing Hamiltonians. This arc is con- 
structed by looking at the cross-sections of H (k) between 
kd = and kd — tt. The rest of the Hamiltonian can 
be reconstructed using inversion symmetry The end- 
points have to be on Id-i because the inversion takes 
the kd = 0,7r cross-sections to themselves. 

Thus, let us solve the following problem: consider arcs 
71,72 in Hd-i connecting two points in the subspace 
Id-i- What conditions ensure that it is possible to move 
arc 71 to arc 72? This deformation is possible if we can 
first slide the end-points of 71 within Id-i onto the end- 
points of 72 and then smoothly deform the curves con- 
necting them. Fig. [3] illustrates the problem. 

We can thus classify d-dimensional Hamiltonians by 
solving two problems: describing the different compo- 
nents of Id—i, and classifying the arcs connecting a pair 
of points in Jid-i up to homotopy. 

Let us now consider d = 3. The first step is analogous 
to the problem we arc trying to solve, just in one dimen- 
sion less. (The components of I 2 are just the different 
classes of 2-dimcnsional inversion-symmetric Hamiltoni- 
ans.) Let us suppose we know the solution to this prob- 
lem, so that the two arcs 71 , 72 can be assumed to have 
the same end-points. 

We now have to slide the interior of arc 1 onto arc 
2. Classifying arcs with fixed end-points in a given 
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FIG. 3: Dimensional induction. The grey region represents 
Hd-i' each point corresponds to a generic d-l-dimensional 
Hamiltonian. The two ellipses on the side represent the com- 
ponents of Td-i, the oHamiltonians with inversion symme- 
try. Inversion-symmetric d-dimensional Hamiltonians (three 
of which are shown) are represented by arcs connecting points 
in Id-i- Two of these Hamiltonians are equivalent if the end- 
points are in the same component of Id-i and have the same 
winding numbers around holes in the space (represented by 
the white ellipse). For example, 72 and 73 are not equivalent 
because their final end-points are in different components; 71 
and 72 are not equivalent because jij^ 1 winds around the 
hole. 

space is closely related to classifying closed loops. The 
loops in %2 can be classified by two winding numbers 
§ da(k z ) and § d/3(k z ), where a and /3 are angular vari- 
ables around holes in "H2- This is a well-known result in 
disguise. Loops correspond to three-dimensional Hamil- 
tonians without any special symmetry, which are classi- 
fied (when there arc sufficiently many bands) by three 
Chcrn numbers. The Chern number of the k z = cross- 
section is a function of the points in not of the loops, 
and it is constant for the whole component we are con- 
sidering, so Ghz does not matter. The remaining two 
Chcrn numbers are the winding numbers. 

Now an arc connecting two points in W2 can wind 
around a hole in the space any number of times just like 
a loop; it just does not close on itself. The topology of 
loops is therefore determined by the value of J Q da{k z ) 
and Jq d(3(k z ) as well. 

The Chern numbers Ghx and Gny of the full Hamil- 
tonian are given by § da{k z ), <f>_^ d(3(k z ); hence the 
"winding numbers" of the open arcs are half as big as 
the Chern numbers, (by inversion symmetry). So if the 
Chcrn numbers of the arcs are equal the Hamiltonians 
are equivalent. 

The only problem left is showing that the end-points 
can be slid to one another under appropriate condi- 
tions. This is the same as classifying inversion symmet- 
ric Hamiltonians in two-dimensions. This problem may 
be reduced to one dimension. The result is that, for two 
Hamiltonians in k x — k y space to be equivalent, the Chern 
numbers Ghz must be the same, and the one-dimensional 
boundary Hamiltonians must be equivalent. 



Now we must classify inversion-symmetric Hamiltoni- 
ans in one dimension. The problem reduces directly to 
zero dimensions because all one dimensional Hamiltoni- 
ans without symmetry are equivalent. 

Two zero-dimensional Hamiltonians (i.e., matrices!) 
are clearly equivalent if the numbers of even and odd oc- 
cupied states are the same-just shift the energy eigenval- 
ues so that the two Hamiltonians match. Hence the last 
condition is that the eight integers n (hL) and the total 
number of occupied states n must match. (The original 
Hamiltonian has bifurcated into eight zero-dimensional 
Hamiltonians through the process of taking boundaries.) 
The number of even and odd unoccupied states above the 
Fermi energy must also be the same, but as mentioned 
above, there are an infinite number of these in continuous 
space. 

Hence, three dimensional Hamiltonians are classified 
by Ghx , Gny and Ghz together with the parities at the 
TRIMs. 



B. Constraint on Parities 

Although the values of the 11 integers n {n) and Ghi 
determine a phase, not all combinations of them are al- 
lowed. First of all, J2 K n °( K ) mus t be even (Eq. (QJ). 
There are two ways to see this. 

For the first explanation, let us understand a more 
general question: consider a Hamiltonian (not necessarily 
a gapped one) that is changing as a function of time. We 
will try to understand what happens when the parities at 
the TRIMs change. The parities at a TRIM change only 
when an even state at the TRIM below the Fermi energy 
and an odd state above the Fermi energy (or vice versa) 
pass through one another. Appendix B shows that, each 
time ti (k) changes by 1 by means of such an interchange, 
a pair of Dirac points appears or disappears. Dirac points 
are defined as points where the valence and conduction 
bands cross, with a cone-shaped dispersion. They are 
stable unless they meet other Dirac points (with opposite 
"chirality") and annihilate. 

Now start from a trivial Hamiltonian, with all electrons 
glued to the Bravais lattice; in this Hamiltonian, all the 
parities are even. After an odd number of changes of 
7t. (k), there are an odd number of pairs of Dirac points, 
so the crystal is not insulating. If the system is insulating 
^2 K n (hi) must be even. 

The stability of the Dirac points is explained in part 
by the basic result on degeneracies of eigenvalues, rather 
than by symmetry: in order to tune a Hamiltonian to a 
point where there is a degeneracy, three parameters are 
sufficient. Since H(k) is a function of three momenta, 
these may be tuned to a point where there is a degener- 
acy, provided H is close enough to having a degeneracy 
in the first place. 

Any set of parities ti {k) satisfying Y[ K Vk = 1 can 
be realized in an insulator. There is never a direct 
phase transition (even with fine tuning) between two such 
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phases when two or more r] K 's flip sign. When two modes 
cross at one TRIM in order to change the value of n 
there, Dirac points will form, and the system will be a 
semimetal. The Dirac points must then move to the sec- 
ond TRIM and reannihilate, so that the system becomes 
an insulator again, as illustrated in Fig. m 

The alternative derivation of the constraint studies the 
Bloch cigenfunctions for a fixed Hamiltonian. Let us first 
suppose there is a single occupied band \tpiw)- To deter- 
mine whether the wave function is even or odd, let us 
take its overlap with an even orbital \s) centered on the 
origin. Define si(k) = (s^ik)- Plot the solutions in the 
Brillouin zone to 

si(k)=0 (14) 

This equation is a complex equation, amounting to two 
equations in three variables, so its solutions are curves. 
At a TRIM, ipi K (r) is either even or odd. If it is odd, 
its overlap with |s) vanishes. Generically, the converse is 
also true. Hence, there is one curve through each TRIM 
where IV'ik) is odd. 

But since the curves are inversion symmetric, they 
must pass through an even number of TRIMs (see Fig. 
[5]). Hence, the total number of TRIMs where Iipi K = 
—ipl K is even. 

When there are several filled bands which do not touch 
each other, Yin 7 ! -^) = 1 f° r eacn band separately, and 
YIk. Vk, = 1 follows. If bands do touch, another step is 
required (the product for a single band may be —1, but 
the total will still be +1). Consider the curves deter- 
mined by s a (k) = (s|V>ak) = for all the occupied bands 
1 < a < n. Some of these curves may be open arcs be- 
cause s a (k) becomes discontinuous at Dirac points. Such 
arcs can only end at a Dirac point where the band touches 
another band. In this case, there is also an arc leaving 
the Dirac point in the other band (see App. [Cj. Putting 
all the arcs from the occupied states together therefore 
produces a set of closed curves; these curves will be var- 
iegated if the arcs in each band are imagined to have 
different colors, but they are still closed. Hence we can 
still deduce that J2 K n o(i^) is even. 

C. Polarization and Hall Conductivity 

The constraint just derived is not the only constraint 
on the 11 invariants. The parities of the Gh arc de- 
termined by the parities of the n 's (see Eq. ©). We 
will prove this momentarily, but it is logically necessary 
to derive the expression for the polarization first. Both 
these quantities are related to the Berry connection, a 
vector function in momentum space. For a single band, 
the Berry connection is defined by 

A a (k) = i(w a k|V k |M ak ) (15) 

and the total Berry connection A(k) is the sum of the 
Berry connections of the occupied bands. 



We will prove that the intrinsic polarization is given 
(modulo e times a lattice vector) by Eqs. (O,©. In 
general, the polarization per unit cell is given by 

P = P e -J2Zie(xi-x ). (16) 

i 

The second term is the polarization of the nuclei in the 
unit cell relative to an origin, xq, and the first term is 
the polarization of the electrons relative to xq ■ Since the 
electrons are delocalized, calculating this contribution is 
subtle. It is given, according to Ref. [lU, by 

Pe = ^E/ dkA ^ k )- ( 17 ) 

n 

This expression for the polarization is ambiguous up 
to multiples of e, as it should be, on account of surface 
charge. If the unit cell is redefined, some nuclei loca- 
tions are shifted by one unit. Likewise if the Bloch wave 
functions are redefined by u n k — > e l8 ^ k 'u n k, then the po- 
larization shifts by j^(9(2ir) — 0(0)), an integer multiple 
of e if e 10 winds around the unit circle. 

To evaluate the polarization of an insulator with inver- 
sion symmetry, set x$ = 0, the inversion center. First, 
consider a single band. The wavefunctions at k and — k 
must be the same up to a phase, so = e^^ijjk 

for some phase 8(k). Therefore, Ak + A-k = 9'(—k). 
Combining k and — k together in Eq. (fl~7|) leads to 

is the parity of the wave function ±1 at TRIMs k = 0,ir. 
Hence if the parities at the TRIMs are different then the 
polarization is a half integer. If there are many bands, 
we may sum the polarization over all of them and we find 
in general that 

(-l) 2p <= m r h . (18) 

The result for the polarization in three dimensions is 
derived from this below, but Fig. QJL suggests the rea- 
son. Each of the four vertical lines through TRIMs looks 
like a one-dimensional insulator with half-integer polar- 
ization, so the net polarization per unit cell of the three- 
dimensional crystal is also IR3. Note that half of any 
Bravais lattice vector is an inversion center, so there are 
eight incquivalent inversion centers. The parities depend 
on which center they are measured relative to, so P e is 
origin-dependent. But P is not because the contribution 
from the nuclei depends on the origin in the same way. 

Now let us consider the Chern number for a two di- 
mensional system, H(k X: k y ) : and show that (— 1) Gh * = 
Y\ K rj K . This is the two-dimensional version of Eq. ([5]) . 
The Hamiltonian leads to a 1-D Hamiltonian when 
k y is fixed. As k y changes, the polarization P(k y ) of 
the one dimensional system changes. This means current 
must flow from one end to the other. According to Thou- 
less's pumping argument, the Hall conductivity Gh z is 
equal to the total charge (divided by e) that flows when 
k y changes by 2tt: (In real space, the one-dimensional 
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FIG. 4: Changing the parities of bands at two TRIMs. The figures represent the Brillouin zone of a system with two bands, one 
of which is filled. Initially, all the filled bands have parity +1, but the parities at two points are changed with the assistance 
Dirac points which also act as monopoles in the Berry flux (see appendix [Cj. A pair of monopoles forms at one TRIM and 
they move to another TRIM where they annihilate. In the process, the parities of the states at both TRIMs are reversed. The 
open circles indicate where the monopoles start out and disappear. 




FIG. 5: Curves which are inversion symmetric must pass 
through an even number of TRIMs. Left, an attempt at draw- 
ing a curve that passes through one TRIM fails to be inversion 
symmetric. Right, an inversion symmetric figure; if there is 
just one curve passing through one of the TRIMs, it must 
go all the way around the Brillouin zone and pass through 
another TRIM on the way. 



system is just the two-dimensional system rolled into a 
tube along the y direction, with period equal to a sin- 
gle cell. Changing k y corresponds to applying an EMF 
around the y-dircction.) 

Thus Grz = —\ Jl 7T dP(k y ). (The polarization is not 

single-valued if G z ^ 0.) Now if Y[ K v( K ) = — 1 ( as m e i~ 
ther of the 2D layers in Fig. [T]b) then the polarizations at 
k y = and k y = ir differ by a half integer. Thus, J* dP = 
(k + i)e. By inversion symmetry, dP(k y ) = dP(—k y ) 
(since P(k y + dk y ) - P(k y ) = -P(-k y -dk y ) + P(-k y )). 
Hence between — tt and the polarization changes by the 
same half-integer, and the full change in the polarization 
is odd, G H z = -(2ft + 1)2tt. 

This section concludes with the generalizations of these 
results to three dimensions. The expression for the Hall 



coefficient in three dimensions Eq. (J5J) is basically a re- 
statement of the two-dimensional result. Each of the 
three components of G# is equal to the two-dimensional 
Hall coefficient for any cross-section of the Brillouin zone: 



Ghz — 



dk z 

(27T) 



G 2d 



Hz 



(kz) 



G 2 Hz(kzo) 



(19) 



since the Chern number for any cross-section k z = k z0 is 
the same since it is a topological invariant. 

The Hall coefficient can be obtained modulo 2 by look- 
ing at an inversion symmetric plane, either k z = or 
k z = tt, giving Eq. ([5]). Note that this gives another rea- 
son for the constraint Y[ K rj R = \: the two planes have to 
agree about G# z 's parity. 

The expression for the Hall coefficient can be guessed 
by looking at Fig. |4] The initial phase is trivial. When 
the monopoles move through the Brillouin zone and an- 
nihilate, they leave behind a magnetic flux, so the Chern 
number is 27r. The flux's direction is parallel to the 
edge where the parities of the occupied states have been 
flipped. 

The polarization in three dimensions is well-defined if 
G# = 0; then there arc no surface modes. According to 
Eq. [3 the pattern of signs is then a plane wave, so all 
the Id polarizations at TRIMs, Pxdi K yi K z), are the same. 
This value in fact coincides with the three dimensional 
polarization. 

The three-dimensional polarization is the inte- 
gral over onc-dimcnsional polarizations, P x = 
N^r P U k v^ z ), if G Hy = G Hz = 0. (This 
condition ensures that P x is single- valued.) The in- 
tegrand is not a constant, but by inversion symmetry 

p? d (k y ,k z ) - Pf d (o,o) = -(P¥ d (-ky,-k z ) ~ Pf d (o,o)). 
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(Intuitively, one expects P* d (k y ,k z ) to be an odd 
function of the wave number, but only differences in 
polarization are inversion symmetric because of the 
ambiguity in the polarization.) Hence P x = Pf d (0,0), 
which is given by Eq. (JT]) and Eq. ©. 

Note that if Gh is nonzero, then part of the polar- 
ization is still well-defined. If Gh is parallel to gi then 
this component of the polarization cannot leak through 
the surface modes. The chiral modes circle around the 
x-axis, so they do not provide a short circuit between 
the right and left faces of the crystal. In fact, the same 
reasoning shows that P x is still given by Eqs. (J?]),©. 

D. Magnetoelectric Response 

Now we justify the relation 9 = ti (k) (Eq. (|6])) 

by using the classification of topological insulators to de- 
form an arbitrary inversion-symmetric Hamiltonian H to 
one that also has time-reversal symmetry. Eq. (|6|) is 
known for Hamiltonians with time-reversal in addition 
to inversion, either by directly evaluating the expression 
for 6 in terms of the Berry connection^ or by combining 
the results of@0and[l 

Hamiltonians with time-reversal symmetry and inver- 
sion symmetry have extra constraints on the Gh and 
n D 's. The states at the TRIMs come in Kramers dou- 
blets, so ti (k) and n are even, and the Hall conductivity 
must be by symmetry. By the classification theory, 
an inversion symmetric insulator can be deformed into 
one that is also time-reversal symmetric if it has these 
properties, too, and Eq. (|6]) follows. 

Now we are already assuming that the Hall conductiv- 
ity is zero, since the gapless surface modes can interfere 
with the magnetoelectric response. When some of the 
n D 's are odd there is an additional step before it is possi- 
ble to deform H into a time-reversal symmetric Hamilto- 
nian. We superimpose an independent inert system onto 
H and then deform the combined system. Consider an 
"ionic" or "frozen" insulator with protons on the Bravais 
lattice and electrons fixed on certain sites, and with all 
hopping amplitudes equal to zero. This insulator does 
not contribute to 9. If the electrons are at half of a Bra- 
vais lattice vector p = ^ an< i its translates, then it is 
inversion symmetric. 

Though this insulator does not contribute to 9, its par- 
ities depend on k. This is enforced by Eq. since the 
insulator has a polarization. Calculating these parities 
directly sheds light on Eq. ([7]). The parities are obtained 
by transforming to a basis of plane wave states, 

ip>k = E e ' !k ' (R+p) i R +p) ( 2 °) 

R 

where |r) represents an electron at an orbital centered 
around r. When k is a TRIM, l\n) p = e 2ip K |/t) p , as- 
suming the orbitals are even under inversion. 

Now consider an arbitrary Hamiltonian % without Hall 
conductivity. Eq. ([7]) implies that the net parities t](k) 



vary just as they do in the ionic crystal (except possibly 
with an overall sign). Therefore, we can add a copy of 
an ionic crystal with a polarization that matches the in- 
trinsic polarization of H, leading to an extended Hamil- 
tonian "Hoxt- All the r^'s will become ±1. This sign 
can be changed to +l-redefine the orbitals of the frozen 
crystal to be odd under inversion. (If all the signs were 
initially —1 add a frozen crystal with p = and odd 
parity orbitals.) 

Now the ti (k)'s have become even. If n is odd then 
add an unpolarizcd crystal, with a single electron in an 
even orbital on the origin and its translates. Now both 
n and u (k) are both even. 

This enlarged crystal's 9 may be determined by de- 
forming it to a time-reversal symmetric insulator, so 
= \ n 0>ex t(/c). But it has the same magnetoelec- 
tric response (9 0X t as the original crystal, and the frozen 
crystals changed the number of odd states at the TRIMs 
by a multiple of 4; hence this formula applies to the orig- 
inal crystal also. 



E. Other Quantized Responses? 

We have found that there are eight integers n (n) that 
are constant within a given phase. However, the calcula- 
tions so far have only found that a few mod. 2 combina- 
tions of them that have interpretations as response func- 
tions. We think that there are no quantized responses 
besides 9 and Gh, for any applied field or geometry. For 
any set of n D (/«)'s, such that Gh and 9 are zero, it is 
possible to find a Hamiltonian where the particles have 
zero hopping. 

Let us start from a general vector of values n {n) and 
try to build it out of the n Q vectors for the frozen (or 
ionic) insulators described in the previous section. Frozen 
insulators do not have any response. Regard n (n) as an 
eight-dimensional vector of integers n Q . The vectors in 
this space can be pictured by labelling the vertices of a 
cube with eight numbers. 

The frozen insulators described in the previous section 
give eight vectors f p in this space, as p varies over all 
possible polarizations, (fo represents the insulator where 
each n = I.) Since n Q is eight-dimensional and there are 
eight fp's, every insulator's vector can be decomposed in 
terms of these frozen vectors. This suggests that every 
insulator can be built out of frozen insulators. This is 
not really true though, because the superposition may 
involve non-integer coefficients. There is no way to take 
a fractional multiple of a system. 

For each vector that can be expanded as n Q = ^ p r p f p 
where the r p 's are integers, there is a frozen insulator 
with the same n Q up to a constant. Hence these n Q 's do 
not require nontrivial quantized responses. For each p =/= 
0, combine together |r p | copies of the insulator with the 
polarization P e = p. (The parities of the inert orbitals 
should be chosen based on the sign of r p .) 

How close can one get to a get n D can one get with 
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such a sum of frozen vectors? Appendix [Gl shows that for 
each eight-component integer vector n Q one can subtract 
a frozen vector f so that 

HQ - f = U X yV X y + U XZ V XZ + Uy Z Vy Z + W X y Z V X y Z (2l) 

where the u's are each or 1 and w xyz is 0, 1,2 or 3. 
The v's are the following vectors. The first, v xy cor- 
responds to the insulator illustrated in Fig. [T|d, with 
the quantum Hall response. It has ones on the edge of 
the cube described by k x = n y — it and zeros every- 
where else. The other u's with two labels are defined 
symmetrically. The vector v xyz has a one at one corner 
of the cube (k x = K y = k z = tt). If the vector n Q of 
an insulator is decomposed in this way, the constraint 
S re n o( K = (mod 2) implies that w xyz is even (either 
or 2 modulo 4). Hence the difference betweena generic 
insulator and a frozen one is characterized by values of 
four numbers-each u can be or 1 and w xyz can be or 
2. Thus any property of the responses that is contained 
in the n D 's is determined by these four parameters. 

Eqs. 1211 and 151 show that the u's describe the quantum 
Hall conductivity modulo 2, while = tt if Gjj = and 
w xyz = 2( mod 4). There seem to be no other responses. 

In spite of this, it is possible to measure some infor- 
mation about f by considering the polarization. (The 
argument does not rule this out because polarization is 
static.) 



IV. PARITIES AND THE ENTANGLEMENT 
SPECTRUM 

The relations between the properties of the insulator 
and the parities can be derived geometrically using the 
entanglement spectrum of the insulator; there is a rule 
for counting the number of states in this spectrum based 
on the parities. 

The entanglement spectrum is defined using the 
Schmidt decomposition. The insulator is cut by an imagi- 
nary plane passing through a center of inversion symme- 
try. The many-body ground state wave function then 
decomposes, 

l*> = 7f E e ~*H;»«- ( 22 ) 

* a 

where Z is a normalization constant and E^ controls the 
weight of a given term. (It is called the entanglement 
"energy" because this formula looks similar to the Boltz- 
mann distribution. The states can be called the entan- 
glement states.) 

When the wave function of the entire system is a Slater 
determinant, the entanglement states \q)l are Slater de- 
terminants, too^. They are formed by selecting wave 
functions from a special family of single-particle wave 
functions /^(r). (These states may be labelled by 
the momentum along the surface, k^, by translational 



symmetry.) Each of these wave- functions has an associ- 
ated "energy" eLi(k±) (as if they were eigenfunctions of 
a single-particle Hamiltonian) . The entanglement "en- 
ergy" E e a is the sum of all the "energies" of the occu- 
pied states. The entanglement spectra, fLi(kj_), can be 
used to determine "topological" properties of the system, 
since these spectra may be continuously deformed into 
the physical surface spectrum of the syste m 11 ! 12 . 

Determining basic properties of the entanglement spec- 
trum is simple in the presence of inversion symmetry. 
The entanglement spectrum has a particle-hole symme- 
try I e that implies a rule for finding the number of entan- 
glement states at each surface TRIM. The I e symmetry 
takes each mode to another mode whose momentum k^ 
and "energy" £i(kj_) have the opposite sign. Let us re- 
gard as the Fermi "energy" ; the state in the Schmidt de- 
composition with the smallest "energy" (i.e., the largest 
weight) is obtained by filling up all states with e < 0. 

At surface TRIMs I e ensures that states appear in 
pairs with energies ±e when e ^ 0. There can be a single 
mode at zero. This mode will stay exactly at zero no 
matter how the system is changed, because moving away 
would break the symmetry More can be said about the 
zero-energy states: the "index" at each TRIM can be 
determined; this is the difference AN e (n±) between the 
number of modes of even and odd parity. This quantity 
is significant because even and odd zero-energy states 
can "cancel" one another and move to nonzero energies 
±e, while two states of the same parity cannot. If two 
states evolve into the eigenstates f\ and fi with energies 
of opposite sign, then f\ and fi are orthogonal states 
exchanged by I e . The corresponding parity eigenstates, 
-^(/! ± / 2 ), have opposite parities. 

The imbalance number can be found directly from the 
bulk band structure, 

AN e (K X ) = 1/2(AAT(ki) + AN(k 2 )) (23) 

where K\ and k 2 are the two bulk inversion-symmetric 
momenta that project to k± and AN(ni) (e.g.) is related 
to n (Ki); it is the difference between the number of even 
and odd occupied states at «i, that is n — 2n {K\). In 
one dimension, for example, there is only one AN e (since 
there is only one surface momentum) to determine, and 
AN e is equal to the number of even states at k = tt minus 
the number of odd states at k = 0. The parities of the 
bulk states are to be calculated relative to an inversion 
center that is on the cutting plane. (Parities depend on 
the inversion center, as does the entanglement spectrum. 
The physical responses described above do nol^ 

The derivation of this requires some results of Refs. 
[20II21I whose derivations are summarized in appendix [D] 
As the structure of the entanglement spectrum suggests, 
the entanglement modes are actually eigenfunctions of a 
Hamiltonian Hl defined on the part of space to the left 
of the cutting plane. The eigenvalues of Hl are not equal 
to 6i, but they are related to them, 

H L \fk ± ) = itanhie L4 (k ± )|/f k J. 
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The Hamiltonian Hl can be obtained from a Hamilto- 
nian in the whole space, Ha at the flat-band Hamiltonian. 
This Hamiltonian has the same eigenfunctions as the true 
Hamiltonian but has different eigenvalues, — A for the oc- 
cupied states and +i for the empty ones. Hl is obtained 
from this by cutting off the right half of the space; cut- 
ting away the left half of the space leads to a partner 
Hamiltonian Hr (whose eigenfunctions fn generate the 
Schmidt states on the right). 

The flat-band Hamiltonian has an unusual property- 
its eigenstates can be reconstructed from the eigenstates 
of the two halves. There is a correspondence M. between 
the eigenfunctions of Hl and Hr which reverses the sign 
of the "energy" . Using this correspondence, define 



MfL(r) ; (24) 



then Fik is an eigenstate of £Zfl a t with eigenvalue —\ 1 
and hence is occupied. As ff 4 varies over all eigenstates 
of Hl, the function F varies over a basis for the occupied 
states in the ground stated 

When the system is inversion symmetric, Ai and T 
can be combined together to give the symmetry I e ; it 
is a transformation within the left half of the insulator, 
defined by I e = lAi. Since 1 is a symmetry of the 
wave function, it preserves e while M. reverses its sign. 
Therefore, I e acts as a particle-hole symmetry 

Now we return to the zero- "energy" states at TRIMs 
and their X e -parities r) e i K± . Each state fi K± extends, by 
Eq. p4|) . to an occupied state 



1 



V2 



(25) 



This state is invariant under I, and the parity is rjein- 

Let us determine the value of AN e for a one- 
dimensional system. The result in higher dimensions fol- 
lows since we can use conservation of kx to reduce the 
dimension to 1. Consider a circular chain with an even 
number of cells, L. Now, count the number of even occu- 
pied states W e minus the number of odd occupied states 
W , using two different bases. W e — W is equal to tr X 
so it is basis independent. 

One orthonormal basis will be obtained by cutting the 
system along a diameter. There will now be two cutting 
points and 4. The zero- "energy" states give parity 
eigenstates centered on each of the two cuts, according 
to Eq. (J35J). These contribute 2AN e to W e — W . The 
remaining states can be organized into inversion-related 
pairs, Fi(x),Fi(—x). All these states are mutually or- 
thonormal because they correspond to different eigen- 
values of Hl,Hr. The inversion matrix I has only off- 
diagonal matrix elements between Fi(±x), so they do not 
contribute to the trace. 

On the other hand, instead of the localized wave- 
functions, we can use the extended Bloch functions, 
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FIG. 6: Determining the Quantum Spin Hall index using the 
entanglement spectrum. The figure compares spectra of a 
nontrivial and a trivial system. The left spectrum for each 
system is the entanglement spectrum, and the right illustrates 
how the surface spectrum might look. In the entanglement 
spectrum, inversion symmetry protects degeneracies at zero 
energy at the TRIMs, allowing one to determine the index. 
But the two sets of spectra can be deformed into one another 
(the difference is probably more drastic than illustrated). Be- 
cause time reversal symmetry produces Kramers degeneracies 
at the TRIMs at all energies, the parity of the number of 
modes crossing the Fermi energy does not change. 



V'a(fcx)- The wave functions at momentum ±k x are ex- 
changed, and the wave functions at the TRIMs contribute 
AN{n) + AiV(O)) to W e - W . Eq. (E2D follows. 

We can now count the edge states of a two-dimensional 
insulator that has both time reversal and inversion sym- 
metry, to obtain the formulae from Ref. [H for the indices 
of topological insulators. We will focus on the two di- 
mensional quantum spin-Hall index, since the three di- 
mensional indices are defined in terms of it. The quantum 
spin-Hall index v is the number (modulo 2) of edge modes 
betweeen and tt. Finding surface states in the entangle- 
ment spectrum is easy because of the extra particle-hole 
symmetry (see Fig. [6]); these states remain when the 
spectrum is deformed into the physical spectrum, by the 
standard arguments. 

Consider the entanglement spectrum created by divid- 
ing the system at y = 0. The index is the parity of the 
number of dispersion curves crossing a line e = const. 
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AN e =0 




AN e =l 



AN e =0 AN e =l 

FIG. 7: A model with a single filled band that has an odd 
quantum Hall effect parallel to R z . a) The parities at the 
TRIMs. b) The entanglement states on the xz-fa.ce of the 
Brillouin zone, determined using AN b (k±) = i(AJV(#ti + 
AN{k 2 )). 



AN e -x ^, e -x 



FIG. 8: Determining the polarization from the entanglement 
spectrum, a) Parities for an insulator with one filled band, 
with a half-filled Fermi sea in the entanglement spectrum. 
b)Possible entanglement Fermi arcs for a cut parallel to the 
xz plane. The arcs surround half the surface Brillouin zone. 
If there are no nuclei on y = 0, this crystal has a half-integer 
polarization in the {/-direction. 



between k x = and n. Consider in particular the axis 
e = 0. Strictly between and tt the axis crosses an even 
number of modes: the crossings come in pairs because 
the spectrum is symmetric under l~I e , which just flips 
the sign of e. (Generically, these states will just mix and 
move off the axis.) 

Therefore only the modes at the ends of the interval are 
important. We may assume that all the modes at either 
of these TRIMs have the same parity, because otherwise 
the states whose parity is in the minority may combine 
with the states in the majority and become gapped. Then 
by Eq. (f2"5| the number of modes at, e.g., k x = is 

|AJV e (0)| = |n-n o (0,0)-n o (0,7r)|. (26) 

Since these modes are at the extremes of the interval from 
to 7r, they only qualify as half-modes. One way to see 
this is to look at a line slightly above the axis. This line 
crosses half of the modes emanating from each TRIM, so 
the number of crossings, mod. 2, is V = 5 X^fe =0 \ n ~ 
n (k x , 0) — n (k x , n)\. This is congruent to \ Y1 K n o(«0, 
summed over all four TRIMs. 

When the flat band Hamiltonian is deformed into the 
true Hamiltonian, v remains the same even though the 
surface states no longer have particle-hole symmetry. 
The energy curves form continuous loops or zigzags (see 
figure) because of the double-degeneracies protected by 
Kramers' theorem. 

We can use a similar approach to understand the re- 
sults for the polarization and Hall coefficient. These ef- 
fects may be determined by sketching the arcs of the en- 
tanglement Fermi surface rather than the real one. While 
each zero energy state at a TRIM gives at least a Fermi 
point, not all of these points extend to Fermi arcs. It 
is certain that the 'parity of the number of Fermi arcs is 
the same as AiV e 's parity however. (This is proved in 
appendix [E]) 

Fig. [7] shows how the modes might look for a set of 
parities that corresponds to an odd quantum Hall con- 
ductance. Consider the xz surface of this insulator. Ac- 
cording to Eq. there must be one arc (or an odd 



number) passing through (71", 0) and (7r, 7r), and an even 
number through the other two TRIMs. Hence as one 
travels along the x-direction through the Brillouin zone, 
one crosses an odd number of Fermi arcs. (Any extra 
arcs that do not pass through TRIMs come in pairs by 
inversion symmetry.) Since the number of arcs (counted 
with a sign depending on the sign of the group velocity) 
is the z-component of the Hall conductivity, G z /(2w) is 
odd. 

Now consider the parities in Fig. [5] One possible 
choice of modes is illustrated. If these modes are not 
chiral, then the states below energy e = cover half the 
Brillouin zone (exactly half because of the symmetry). 
This will imply that P y = | if there are no nuclei on 

= 0. 

The polarization is defined as the surface charge once 
the surface bands have been emptied. (See appendix [F]) 
This definition may be applied to the entanglement spec- 
trum. Cut the system with a plane y = through a point 
of inversion symmetry. An important property of the en- 
tanglement spectrum is that every term in the Schmidt 
decomposition has a net charge of zero in the vicinity of 
the cutting surface. 

Of all the terms in the Schmidt decomposition, con- 
sider the one with the largest coefficient, \G)l\G)r, 
where \G)l and \G)r are the ground states of Hl and 
Hr respectively, obtained by filling all the negative en- 
ergy states. (It does not matter which of the states with 
exactly zero energy are filled.) These states are mirror 
images of each other, so they both have the same surface 
charge Q y . The net charge is 2Q y . However, the Schmidt 
state \G)l\G)r, obtained by a hypothetical measurement 
of the two halves of the system, has the same charge as 
the ground state. Each electron is pushed over to the 
side of the y = plane where it is more likely to be, but 
no electrons appear or disappear. Thus, Q y = 0. 

Now the ground state of Hl has a partially filled band. 
When this band (which covers half the Brillouin zone) 
is emptied, a charge of — | per unit cell on the surface 
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remains. Hence the polarization P v = | (mod e). 

To compare this result to Eq. ([8]), we need to take 
nuclei into account. Note that Pf = for the given 
parities. By neutrality, there must be a single charge 
— e nucleus per unit cell, which contributes P% = —%. 
Appendix IF1 gives more details. 

Conclusions This study of inversion symmetric insula- 
tors displays various interconnections between quantized 
electrical properties and the inversion parities at TRIMs. 
In particular, it provides a simple criterion for determin- 
ing when a crystalline compound has a magnetoelectric 
response equal to it. The magnetoelectric effect could 
be seen in a system with broken time-reversal symmetry, 
even if the spin-orbit coupling is small. 

The properties of inversion symmetric insulators can 
be derived in a simple way with the help of the entangle- 
ment spectrum. 

Future work could extend these results to other inter- 
esting lattice symmetries, as well as the classification of 
superconductors with crystal symmetries. Finally, un- 
derstanding the stability of the parity invariants to in- 
teractions will be an interesting subject for future study. 
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Appendix A: Classification using the topology of 
Hamiltonians 

In this section we will show how to use systematic 
methods from topology to classify the insulators with in- 
version symmetry This method was also used to find the 
phases of systems with one of the Altland and Zirnbauer 
symmetry groups^. 

In general we write a band insulator as -ff(k), which 
we think of as a vector bundle on the Brillouin zone. 
The eigenvectors of H(k) with energies below the chem- 
ical potential are the occupied states, defining a vector 
subspace of the entire Hilbert space at k. For inversion- 
symmetric insuators, there is an additional constraint on 
the Hamiltonian: P°[H(k)} = I H(k)I = H(-k), this 
relates the Hamiltonians at k and k. 

The idea behind our method of classification, similar 
to the one used for time-reversal invariant topological 
insulators, is to only look at half the Brillouin zone. 

For a Hamiltonian H in e?-dimcnsions, one can con- 
struct a path f(t) = H(k x = 7rf, k 2 , ^3, • ■ • ) in the space 
of (d— l)-dimensional Hamiltonians, where the endpoints 
(t = 0, 1) are inversion symmetric. Thus the classification 
of d-dimensional Hamiltonians is equivalent to the classi- 
fication of paths in ^(Hd-uld-iild-i), where Hd is the 
set of general d-dimcnsional Hamiltonians and Id C "Hd 



is the subset which is inversion-symmetric. The symbol 
<P(A; A, B) is defined to be: 

«P(A; A, B) = Set of paths in X from A to B (Al) 
= {f:[0,l]^X\f(0)eAJ(l)eB}, 

where A, B C X. We want to divide up Id into sets such 
that paths from different sets are not nomotopic to each 
other. Heuristically, the classes of paths of ty(X;A, B) 
are given by the number of components of A, B which 
determines the set of possible endpoints, and the loop 
sturcture of X which determines the number of ways to 
travel from A to B. The classes of insulators in d + 1- 
dimensions is very roughly given by: 

Components of Id+i 

Loops in U d , n > 2 ... 

~ x (Components olid) ■ (A2) 

Loops in Id 

This idea is made precise using algebraic topology^, and 
is captured by the exact sequence (|A4[) . 

There is an addition structure in the classification of 
insulators^, which comes from the fact that one can com- 
bine two insulators together using direct sums "©" . To 
simplify the classification, it is useful to also have a sub- 
traction "0" operation between insulators. This would 
give the topological invariants {e.g. Z) a group structure. 

The subtraction procedure is realized by considering 
an ordered pair of bands (H\,Hii), which represents 
the "difference" of the two Hamiltonians. Addition by 
H' is given by (Hi © H',H 2 ) and subtraction is given 
by (Hi,H 2 © H'). Imposing the equivalence relation 
(Hi © H', H 2 © H') ~ (Hi,H 2 ) makes the addition and 
subtraction processes cancel each other. Physically, we 
are interested in classifying difference of two topological 
insulators - this is analogous to studying domain walls 
whose properties are determined only by the difference 
in topological invariants. With this interpretation, it is 
possible to talk about a negative number of filled bands 
(whenever H 2 has more bands than Hi). 

The construction above, called the Grothendieck group, 
yields two direct results. First, two insulators Hi,H 2 are 
deformable to one another when the toplogical invariants 
of the band structure (Hi,H 2 ) are all trivial. Second, 
when the latter insulator is the vacuum, the topological 
invariants of (H, vac) are simply the bulk band invariants 
of H. 



1. Hamiltonians, Classifying Spaces and Homotopy 
Groups 

Consider an x TV matrix H with n occupied states 
and N — n empty states. Setting the chemical potential 
to be zero, H is a matrix with n negative eigenvalues and 
N — n positive eigenvalues. In the topological classifica- 
tion of insulators, the energies are irrelevant so long as we 
can distinguish between occupied and unoccupied states, 



and hence we can deform the energies (eigenvalues) of 
all valence bands to —1 and the energies of conduction 
bands to +1. Wc can also assume there are an infinite 
number of conduction bands, and so we let N — > oo. Ho 
is the set of such OD Hamiltonians, and can be separated 
into discrete components based on the number of filled 
bands (which may be negative). The space Ho is home- 
omorphic to Z x BU, where BU is the classifying space 
of the unitary groups. 

At the TRIM, the Hamiltonian is inversion-symmetric 
and commutes with the operator I . The Hilbert space 
is divided into an even subspace and an odd subspace, 
based on the inversion eigenvalues. Hence the set of 
inversion-symmetric (OD) Hamiltonians Iq is homcomor- 
phic to Ho x "Ho- 

Finally, we introduce the "vacuum" vq G Iq C Ho, 
which is a Hamiltonian with no filled bands, vq is a useful 
object in that it allows us to compare any Hamiltonian 
to it, and also acts as the basepoint when we compute 
the homotopy groups of Hq,Iq. 

Given a topological space X and a basepoint within 
the space xo, the homotopy group tt s (X) is the set of 
equivalence classes of maps / : (S s , bo) — > {X, x ), where 
the basepoint b G S s and /(£>o) = x o- For example, 
tto(X) gives the number of connected components of X, 
tti(X) tells us which loops in X are equivalent and which 
loops are contractible. 

The group strucutre of n s (X) is given by concatena- 
tion of maps. However, the group structure of the Ham- 
litonians has already been defined based on direct sums. 
Fortunately, the group composition defined based on the 
two methods (concatenation / direct sums) are compat- 
ible. 

The homotopy groups of Hd are known: 7ro (Ho ) = Z 
because O-dimensional Hamiltonians are classified by the 
number of filled bands n. tt\(Ho) = means that the 
loops in Ho are all contractible. ^("Ho) = Z, because 
maps of the sphere arc classified by the first Chern class 
(or the Chern number) . This invariant give rise to the in- 
teger quantum Hall effect. For higher dimension, -k s (Ho) 
is when s is an odd integer, and Z when s is an even 
integer, corresponding to the (s/2) th Chern class. In this 
section, Chern numbers are taken to be integers rather 
than multiples of 2tt. 

The homotopy groups of To are simply the squares 
of the homotopy groups of Ho- In particular, the set 
of components 7ro(Zo) = T? is labeled by two integers: 
(n, a), where n is the total number of valence "bands" 
and a = n Q is the number of states which have odd in- 
version parity. 



2. Relative Homotopy Groups and Exact 
Sequences 

The homotopy groups n s (X) classifies components, 
loops, and maps from higher dimensional spheres to the 
space X. The relative homotopy groups tt s (X, A) classify 
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FIG. 9: Relative homotopy groups. The figure represents an 
element of TV2(X, A, xo) where X is R 3 ; A is the torus. 

maps from paths, disks, etc. where the boundary must 
lie in some subspace (see Fig. This is how topolo- 
gists define "winding numbers" of open arcs, which were 
discussed in Sec. IIII A[ Relative homotopy groups were 
applied much earlier by Refs^ to an interesting prob- 
lem within physics: classifying defects of ordered phases 
when the defects are stuck to the surface. 

Given a space X a subspace A C X, and a basepoint 
.t G A, the relative homotopy group ir s (X,A) is the 
equivalence classes of maps (D s ,dD s ,bo) — > (X,A,xq). 
The boundary of the disk dD s = must map to A, 

and the basepoint 6o 6 dD s maps to xq. The relative ho- 
motopy groups can be computed via the exact sequence 

TT S {A) ^ 7T S (A) ±> 7T S (X,A) A 7T,_l(A) ±> 1T s —l(X). 

(A3) 

An exact sequence is a sequence of groups along with 
maps defined from one group to the next, each map pre- 
serving the group operations. In an exact sequence, the 
image of every map is the same as the kernel of the sub- 
sequent map. 

In the one-dimensional case, the relative homotopy 
group tti (X, A) describes the set of possible paths in X 
from xo to A up to homotopy, that is, the classes of paths 
*P(A; xo, A). The exact sequence becomes: 

tti(A) ^> Tripo ^> MX, A) A ma) ^> MX). 

(A4) 

In the exact sequence above, the maps are defined as 
follows. 

• i : A — > X is the inclusion map which takes every 
point from A to a point in X. The induced maps 
i( s ) : ir s (A) — > 7r s (A) and takes components/loops 
from one space to the other. 

• All the loops in X start and end at xo, and so they 
are also clearly paths in *J?(A; xo, A) seeing xo £ A. 
jn\ is the map that takes the equivalence classes 
of loops tti(X) to the equivalence classes of paths 
MX, A). 

• d : 7Ti (X, A) — > 7To (A) is a map that takes a path 
and selects its endpoint to give a component of A. d 
is called the boundary map, it takes a "ID object" 
to give a "OD object." 
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FIG. 10: Understanding the exact sequence. Paths in 
7Ti (X, A) can be enumerated by taking one path for each 
inequivalent end-point and adding loops from -k\(X) to the 
paths. This description is sometimes redundant. In this fig- 
ure, X is the cylinder and A is the two circles at its ends. A 
path p and a path p' = p + l obtained from it are shown where 
/ winds around the cylinder. The two paths are equivalent be- 
cause they can be deformed into one another by bringing the 
right end-point around A. This happens because I can be 
smooshed into A. 

It appears that the maps i and j "do nothing" to the 
objects (points, loops) they act on. However, each map 
gives the loop/path more freedom to move around. For 
example, jry takes a loop to a path where the endpoints 
no longer have to be the same. 

The exact sequence captures the idea that the paths in 
ty(X; xq, A) can be classified once one knows the proper- 
ties of X and A, based on their end-points and how they 
wind. 

1. First, we pick the endpoint x G A of the path p. x 
must be connected to xq- and this is captured by 
the statement ker(i( O = img(<9). The explanation 
is as follows. The endpoint of the path is given by 
the d map: x = dp. The image of d is thus the 
set of elements of ito(A) that are connected to .t 
within X. On the other hand, these are exactly 
the components which become equal to in tto(A'), 
when A is enlarged. Hence ker(i( )) = img(<9). 

2. Given a choice of a path p from xq to x, we can 
construct all the other paths between the points. 
We can create any other path p' by concatenating 
a loop I G Ki(X) at the beginning of p. 

3. However, the paths p and p' are only different (i.e. 
not homotopic to each other) only if the loop / can- 
not be unwound within A, this is to say that p ~ p' 
if I is homotopic to a loop in A. See Fig. [TUJ This 
idea is captured by the exactness at 7Ti(A'). Hence 
we think of I belonging to the quotient , 
and this group is called the cokernel of the map 

We see that any path may be constructed by its end- 
point x and a loop I. 

p = j(l) + d~ 1 x. (A5) 

The inverse boundary operator d^ 1 is not unique, but 
that does not affect the structure of the group 7Ti (X, A) 
for the cases we are considering. The equivalence classes 



of x form ker(i( )), while the equivalence classes of I make 
up coker(i(!)). The relative homotopy group is a scmidi- 
rect product 

tti(X, A) = coker(i(x)) x ker(i( )). (A6) 

What this means is that coker(i( 1 )) is a normal subgroup 
of tti(X, A) and ker(i( )) is the quotient of the two. For 
the purpose of classifying inversion-symmetric insulators, 
we can treat the semidirect product as simply a product. 



3. One Dimension 

In this section we examine the classification of ID 
inversion-symmetric Hamiltonians H(k). Let X\ be the 
set of maps H : S 1 -> H Q such that V°[H(k)} = H(-k). 
X\ is homcomorphic to the set of paths in Jio that start 
and end in Xo (i.e. X\ w *P(Ko;2o,2o)), and we seek to 
classify all such paths - to compute 7To(2i). 

For a ID band structure H(k), we can 'factor' out H(0) 
by letting H'(k) = H(k) 9 H{0) such that H'(0) = v . 
H(0) is an element of Ho while H'(k) is an element of 
X\ , where X\ is the subset of X\ with a fixed basepoint 
(iJ'(O) = d ). The decomposition H(k) = H(0) © H'(k) 
can be expressed as 

X 1 =X xl 1 . (A7) 

Hence the classification of X\ can be broken up in to two 
parts, the classification of Xq and that of X\ . The former 
is understood as 7To(2o) = 

Notice that X\ is homcomorphic with the class of paths 
?P(Ho',vo,Xo), whose components arc described by the 
relative homotopy group t:i(%q,Xq). These correspond 
to the possible phases of band structures in X\. The 
relative homotopy group can be computed by the exact 
sequence (|A4[) . Using the fact that 7r (Xi) = ■k\(^H.§,Xq), 

(A8) 

Since tti(Ho) = 0, we can ignore the left side of the 
exact sequence (cokernel of is zero). The integer 
n G tto('Ho) = Z tells us how many filled bands there 
are, and (n , n e ) G tvq{Xq} = Z are the number of even- 
parity and odd-parity bands. The map Z(q) is given by 
n = n + n e , and so the kernel of the map is the sub- 
set where n = —n e , isomorphic to Zi21 It is clear that 
7To(Zi) is isomorphic to ker(i( )) = Z. This is to say that 
the set of paths ^(Ho;vo,Xq) are solely classified by the 
endpoint. Since X\ = Xq x X\. we have 

ir (Xi) = n (X ) x 7T (Ii) 

= 7T (2o) x 7TiC%o,Zo) 

= Z 2 x Z. (A9) 
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The invariant 7To(Zo) = ^ 2 corresponds to the number 
of total bands and odd-parity states at k = 0: n, n o (0). 
The invariant ~kq{Iq) = Z corresponds to the difference in 
number of odd bands at k — it and k — 0: a x — n (jr) — 
n o (0). Hence the parities at the two TRIMs completely 
classify all ID inversion-symmetric Hamiltonians. 

A generator of a group is an element which gives the 
entire group by group addition and subtraction, for ex- 
ample 1 is a generator of Z. In our case, the generators 
arc Hamiltonians which forms a basis for all Hamiltoni- 
ans, up to homotopy. 

The generators of 7To(2o) = Z 2 are two Hamiltonians 
H® and H® which adds one even-parity and odd-parity 
band to the system, respectively. Explicitly: 



so that 



K =[-!](_) QH° n - 



(A10) 
(All) 



where the subscript (±) labels the inversion eigenvalue(s) 
of the orbital(s). The first expression adds an in- 
ert band to increase n; the second expression adds 
an odd-parity band but subtracts an even-parity one 
to increase n Q while maintaining n. The generator for 
7To(2i) = Z is: 



H'Jk) 



- cos k sin k 
sin k cos k 



(A12) 



(+- 



The subscript (-1 — ) specifies the inversion operator Iq = 
[ +1 _J for this Hamiltonian. When k = 0, the matrix 

becomes [ _1 J ^ n and the occupied band is even under 

inversion. Similarly, the matrix at k = it is ) 

and there the occupied band is odd. Hence a x = n (jr) — 
fto(0) = 1 — = 1 and H\(k) is a generator of no(Xi). 

Therefore, any ID inversion-symmetric Hamiltonian is 
nomotopic to 

H{k) = nH° © ctHl © a x Hl{k). (A13) 



4. Two Dimensions 

We apply the same ideas used in ID to classify 
inversion-symmetric insulators in 2D. The inversion- 
symmetric Hamiltonians {k x ,k y ) in 2D satisfy: 
V°H(k x ,k y ) = H{-k x ,-k v ). The set of 2D inver- 
sion symmetric Hamiltonians (Z 2 ) is equivalent to 
*P("Hi;2i,2i), where "Hi is the set of ID band structures 
(loops in "Ho). 

Just as in the ID case where we decompose H (k) into a 
0D and ID object: H(k) = H(0) © H'(k), we decompose 
the 2D Hamiltonian into 0D, ID and 2D components. 
Let 



H'(k x ,k y ) = H(k x ,k y ) e H(0,0), 
so that H'(0,0) = vq. Now we define 



(Al4a) 



H"(k X , ky) = H'{k X ,ky) © H'(0, ky) Q H' \k X , 0), 

(A14b) 



H"{0,ky) =H"(k x ,0) =v . 



(A15) 



H(0, 0) is an element of 1$, and H'(0,k y ) and H'(k x ,0) 
are elements of T\ . We define I2 to be the set of inversion- 
symmetric Hamiltonians satisfying (|A15j) . With this pro- 
cedure, we have decomposed I2 as 



I 2 =I(jxlJx X 2 , 



(A16) 



where Xq = Iq. Explicitly, the decomposition is: 
H[k x , ky) 

= ff (0, 0) © H'(0, k y ) © H\k x ,0) ® H"{k x , k v ) (A17) 



II Xli 



T 2 



Due to (|A15j) . we can think of Hamiltonians in I2 as 
maps from the sphere (rather than the torus) to Jio- 

We now analyze the properties of Hamiltonians in 
X2. For each fixed k x , the Hamiltonian H"(k y )\k T is 
a map from the ID Brillouin zone (S 1 ) to Tio (where 
we've also defined Ho = Ho)- Denote the set of such 
maps as Hi, the set of loops in Ho with basepoint ^o: 
Hi = ty(Ha;vo,vo). Hi is called the loop space of Hi, 
and the notation used in literature is Hi — flHo- 

The Hamiltonian at k x = it is inversion-symmetric, 
and so H" (fcy)^^ 6 Ij. At k x = 0, the line 
H"(ky)\k x — vo is a constant map - which we call vi 
(a line of vo). Clearly vi is an element of Xi, and acts as 
the basepoint when we compute the homotopy groups of 
Xi,Hi- 

Having defined the spaces l2,Hi and basepoint v\, we 
can see that I2 is homeomorphic to the set of paths in 
Hi with endpoints at vi andli: I2 ~ ^P(Hi; vi,Ii). The 
exact sequence which gives the equivalence classes of such 
paths is 



7Ti(Ii) ^ TTi(Hl) ^ 7ro(22) 7T (Xl) 



7ro(«0. 
(A18) 



Elements of tti(Hi) are two-dimensional Hamiltonians 
that equal vq along k x = ~ 2ir and k y = ~ 2ir, which 
are equivalent to maps S 2 — > Hq. Hence 

tti(Hi) = 7r 2 (Ho) = Z (Chern number). (A19) 

The map j , which essentially maps general two- 
dimensional Hamiltonians to inversion symmetric ones, 
is defined by 

(j'HXk^ky) 



< k x < 7T 



V°H{Att - 2k x , 2ir - k y ) ir < k x < 2ir 



(A20) 



It builds an inversion symmetric Hamiltonian out of two 
copies of a Hamiltonian with no special symmetries. The 
map d : I2 is defined by 



[dH]{k)=H{k X ^TT,ky = k) 



(A21) 
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FIG. 11: Isomorphism between ft%(j4a,%s) and 7Ti(Xi). An 
element in Ii is a path from vq to an element of To (dark black 
lines). An element of n 1 (Ti) is a family of such paths. The 
first path in the family must be vi, a "zero- length" constant 
path at uo (gray dot). The last path of the family must also 
be lu, and so the paths must shrink back to vo at the end, 
because that is the base point for the homotopy groups. This 
family of paths traces out a disk D 2 which maps to Ho, the 
endpoints of each path traces out a circle S 1 = 3D 2 which 
maps to lo (dashed circle) . The basepoint of the circle maps 
to Wo, and hence every element of 7Ti(Xi) is also an element of 
7T2(Ho,2o) and vice- versa. This establishes the isomorphism 

7T 2 (Ho,2o) = 71"! (II). 



which takes a 2D Hamiltonian and picks out the ID 
Hamiltonian at k x = n. 

To solve the exact sequence earlier (|A18|) . we need to 
compute 7Ti(Ii), which classifies loops of ID inversion- 
symmetric insulators (not to be confused with 2D insu- 
lators). As argued in Fig. 1111 the group 7Ti(Ii) is isomor- 
phic to 7T2(Ho,2<))- We can compute 7Ti(2i) = ^(HqjIq) 
via the exact sequence (|A3|) . 



z 2 ~* z -^Mii)^ Q ~> o 

(A22) 

The first map i%\ : 7r2(Io) ~~ ^ ^2(^0) is given by (3 = 
fi e + p o (which means that the Chcrn numbers of the odd 
and even bands add to the total.) The image of is all 
of i® 2 y so its cokernel is trivial. Since on the right side 
the groups are also trivial (ker(i° 1 ^) = 0), we have that 

7T 1 (Xi) = 0. 

We now return to the exact sequence (| A18|) to compute 
the tt (1 2 ). 

7Tl(±l) < 7T (?i 2 ) il (f v a 7To(Xi) 7r (Kl) 

Z 01 2J Z ' 

(A23) 

The exact sequence yields 710(22) = Z x Z; the set of 
H"(k) are classified by two integers {a xy ,/3 xy ). The first 
integer a xy coming from the map d gives n (7r,7r); the 
second integer j3 xy is related to the Chern number of H". 
The basis of ^0(12) are found by taking the image of 
the generator in 7Ti(Hi) and one of the preimages of the 



generator element in wq(Ii), which we denote by -ffj(k) 

and iJ 2 (k) respectively. Any Hamiltonian in I2, up to a 
homotopy, can be written as 



H" = (3 xy H 



a xy H a , 



(A24) 



The generator _ff|(k) = i /-n -ff chern (k) where i? c horn is the 

generator of tti(Hi) = ^{Ho), a 2D band insula- 
tor with Chern number +1. H 2 (k) is defined such that 
(dH 2 )(k) = H^{n,k) is the ID insulator ()AT2l) . 
Explicitly, the generator 



1 



1 + x 2 

OK, 

where x = cot % 



1 — x 2 — y 2 2{y + ix) 
2(y — ix) x 2 + y 2 - 1 



(+-) 
(A25) 



1 y 



cot -j-. When x or y = 0, then 
x 2 + y 2 — > 00 and £f 2 (k) = [ _1 J (+_) and n (k) = 0. 

When x = y = n, it can be seen that -ff 2 (k) = f 1 _J ^ j 

and the filled band is odd under inversion. Note that 
the Hamiltonian has Chcrn number +1. Since a xy = 
n (7r, 7r) (and the other TRIMs have n Q = because of 
the normalization), this is related to the constraint G = 
E« "«•(«)■ 



The other generator H 2 (k) 



'(i) 



^chcrn(k): 



Ha{k x , ky) — H a (2k Xl ky). 



(A26) 



has Chern number +1 in each of the halves < k x < tt 
and 7T < k x < 2tt. Evidently, the Chern number of the 
entire Brillouin zone is given by 



G 



2ft X y 



(A27) 



since H 2 , H 2 has Chcrn number 1 , 2 respectively. 

The decomposition (|A16|) gives us six Z invariants in 
2D: two from Tq, one from each of the two Zi, and two 
more from T2- The six invariants (n, a, a x , a y , a xy , (3 xy ) 
are related to the properties of the original Hamiltonian 
H(k x ,k y ) as follows: 

• n gives the number of filled bands, generated by 
(lATUl) . 

• a = n o (0,0) is the number of odd-parity states at 
(k x ,k y ) = (0,0), generated by (jAllj) . 

• a x = ti (tt, 0) — n o (0, 0) involves the difference be- 
tween two parities, generated by (|A12[) (k — > k x ). 

• a,y = n o (0, 7r) — n o (0, 0) involves the difference be- 
tween two parities, generated by (|A12|) (k — > k y ). 

• a xy = n (ir, n) - n (ir, 0) - n o (0, tt) + n o (0, 0) in- 
volves the parities at all TRIMs, generated by 
(|A~25]) . 

• ftxy relates to the Chern number: G = 2j3 xy + a xy , 
generated by (|A26|) . 

The rule for the Chern number's parity, Eq. fol- 
lows from the last constraint here. 
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5. Going to Higher Dimensions 

In d-dimensions, we want to calculate 7ro(X<j), the set of 
components of the space of inversion-symmetric Hamil- 
tonians. To calculate this in a larger dimension, we will 
need to know 7T s (Xd- s ) in lower dimensions. 

The relevant spaces are T-Ld and Id- The general Hamil- 
tonian space Hd = ClHd-i is the space of loops of Hamil- 
tonians in Hd-i, which become trivial on k x = as well 
as on all the other boundaries ki = of the Brillouin zone. 
Their homotopy groups arc given by iT s {%d) = 7r s+( i('Ho). 

The homotopy groups of the subspace Id C Tid are 
harder to find. This space is homeomorphic to Id ~ 
^("Hd-i; Vd-i,Id-x)i since half of the Brillouin zone de- 
termines the Hamiltonian. The homotopy groups of 
Id are therefore given by the relative homotopy groups 
ir s (Id) ~ TT s+ i(Hd~i,Id~i)- Via the relative homotopy 
exact sequence (|A3[) . one can relate 7r a +i['H < i—i,Id—i) to 
7r s ,7r s+ i(%_i) and7r s ,7r s+ i(I d _i). The homotopy struc- 
ture of Id depends on that of Id-i- Iterating this process 
reduces the problem to that of the basic objects Ho and 
X . 

Specializing to 3D, we follow the same prescription as 
before to decompose Hamiltonians. 



x If x X| x I- 



The homotopy group 71-0(13) 
puted, 



(A28) 

7Ti (7^2,^2) can be com- 



7Tl(X 2 ) 



7Ti ("H 2 ) A 





7To(X 3 ) 



d 7ro(22) 2*. 7T (H2) 



z- 



z 



(A29) 



We know that 7To(X 2 ) = Z 2 from the previous section and 
7roC% 2 ) = n 2CHo) = Z corresponds to the Chern number 
G. In addition, the homotopy group iT2{'H.2) = tt^T-Lo) = 
is trivial, so 7Ti(X2) is irrelevant to the problem. 

The map i 2 ^ : 7r (X 2 ) -> ^0(^.2) is given by (|A27[) 
which is surjective, hence coker(i 2 ^) = Z. It follows that 

71-0(23) = Z; 3D insulators are described by the integer 
&xyz- To relate a xyz to bulk quantities, we examine the 
maps d and i 2 under further scrutiny. The relevant part 
of the exact sequence is: 







tto(X 3 ) 9 
Z[a xyz ] Z 2 



7ro(X 2 



'tr- 



J Pyz 



(o). 7ro("H2) 
Z\G] ' 



(A30) 



where the map i 2 ^ is given by G 



-2(3 yz . The kernel 



of the map is the set (a yz ,f3yz) = (2j, —j) for integers j. 
Since the kernel is isomorphic to 770(13), we can identify 
j as a xyz \ hence d is defined by da xyz = {ot yz ,P yz ) = 

{2(X X y Z , &xyz)- 

In terms of the band structure invariants, we have 



2a 



xyz 



n (7r,7r, 7r); 



note that all other n 's are because the Hamiltonian is 
normalized. Explicitly, the generator is 



Hl= [t Q T*+T x {t-<r) 



(++—) 



e 2H2, 



(A31) 



where 



(^0 itxityi^z) — 

with X 



(1 



2y, 2z, 2x) 



cot^Sy 



+ y + z 

cot -^,z = cot %. 



At k x = k y = k z = 7T, (t ,t x ,t y ,t z ) = (1,0,0,0) and 
there are two filled bands with odd parity. At the plane 
k x = tt the Hamiltonian reduces to two copies of (|A25|) . 
but with opposite Chern numbers so that the net Chern 
number is 0. One can think of the two subspaces of dH^ 
as Hi (Chern +1) and /X 2 e Hj (Chern -1). 

Therefore, there are 12 Z invariants in 3D. They relate 
to the bulk qualities as follows. 

• n is the number of filled bands. 



a = 


n o (0,0,0). 




OLx ; 


= n o (x,0,0)-n o (0,0,0), 




a y -- 


= n o (0,7r,0)-n o (0,0,0), 




a z -- 


= n o (0,0,7r)-n o (0,0,0). 




Otyz 


= n o (0, 7T, 7r) - a v - a z - 


0, 


a Z x 


= n (7T, 0,7r) - a z - a x - 


0, 


&xy 


= n (n, 7T, 0) - a x - a y - 


a. 


Pyz 


— 2"( ^ ^2/2)5 




Pzx 


= \{G zx -a zx ), 




fixy 






Oi X y. 


z = 5(n (7T,7r, 7r) - Y^ a ^ 





The last equation explicitly written out as a function of 
n (tt) is 

2a xyz = n (7T,7r, 7r) 

- n o (0, tt, tt) - n (7T, 0, tt) - n (iT, tt, 0) 
+ n (7r, 0, 0) + n o (0, tt, 0) + n o (0, 0, tt) 
-n o (0,0,0). 

From the formula, it is clear that the sum of parities of 
n («) at the eight TRIM is even, Eq. dJ). 

In every higher dimension d, Id has a Z invarant corre- 
sponding to the inversion parity n (rr, . . . , tt) generated 
by H^. For the even dimensions d = 2s, there is a second 
Z invariant corresponding to the s th Chern class, gener- 
ated by Hp. Hence 7ro(X 2s ) = Z 2 and 7ro(X 2s +i) = Z. 

The generator for the s th Chern class is as follows. Let 
{r 1 , T 2 , . . . , r 2s+1 } be 2 s x 2 s gamma matrices satisfying 
the Clifford algebra PP + FT* = 2<5«. First we warp 
the Brillouin zone to a sphere: T 2s [ki, . . . , k 2s ] —> S 2s [h] 
in an inversion-symmetric way by sending k = (tt, . . . , tt) 
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to h= (1,0,0, 
louin zone (ki 
defined as 



. ) and all the planes bounding the Bril- 
0) to (-1, 0,0,...)- The Hamiltonian is 



Hf(k) = ri-r, 



(A32) 



where T is the (2s + l)-vector of gamma matrices. (|A25[) 
is an example of this construction for d = 2s = 2. The 
inversion matrix is given by the first gamma matrix: 
Io = r 1 , and we can see that there are 2 S_1 odd-parity 
occupied states at k = (it, . . . , it). At the other TRIMs, 
the occupied states have even inversion-parity. 

The s th Chern number C s may be computed by the 
formula: 



Tr[P(dP) 



2s 1 



(A33) 



where -P(k) = |(1 — H(k)) is the projector onto the filled 
states and d is the exterior derivative in the Brillouin 
zone. Evaluating the integral shows that C s — ±1 for 
the Hamiltonian H 2s . 

In 2s dimensions, Eq. (|A23|) . generalized to more di- 
mensions, gives a preliminary way of choosing the gener- 
ators: Hp s is the image under j of a generator of itQ(fi. 2s ) 
and H 2s is an arbitrary preimage under d~ l of the gen- 
erator of 7i"o(X)2s— l- Any Hamiltonian can be expanded 
as 



H — fi 2s H 2 ? 



<a 2s Hl s . 



(A34) 



H 2s can be used for the generator H 2s . To see this, we 
decompose H 2s in terms of the original pair of generators 
Hl s ,Hj s . The s th Chern number of Hf is 2, since it 
is constructed using the map j which takes a general 
insulator to an inversion-symmetric one by duplicating it 
in each half of the Brillouin zone [see Eq. (|A20[) ]. At all 
the TRIMs, H 2s (k) = vq and so n a is zero. Using Eq. 
(TA34)) for H 2s implies 

C S {H 2 C S ) = 1 = 2(3 2s + a 2s C s (H 2 a s ), (A35a) 
n (H 2s ) = 2 s - 1 = + a 2s n {H 2 a s ). (A35b) 

The first expression requires a 2s to be odd, and the sec- 
ond requires it to be a factor of 2 S_1 . Hence a 2s = ±1 
and we can use H 2s as the generator H 2s . Since every 
Hamiltonian can be expressed by Eq. (|A34[) . the num- 
ber of odd inversion-parity states n (7r,7r, . . . ) is always 
a multiple of 2 S_1 . 

In terms of the general Hamiltonians in I 2s , the total 
number of odd parity states at all the TRIMs must be a 
multiple of 2 S ~ 1 . Furthermore, 



n (h) = C s (mod 2) in 2s-dimensions. 



TRIM k 



(A36) 



In 2s + 1-dimensions, the Hamiltonians at k x = and 7r 
are both 2s-dimensional inversion-symmetric Hamiltoni- 
ans, and they must have the same Chern number, so the 



constraint on the parities is 

u (k) = (mod 2 s ) in (2s + l)-dimensions. 

(A37) 



TRIM re 



This is related to the (2s + l)-dimensional Chern-Simons 
integral 6 2s +i £ [0, 2tt) by 

1 n (K) (mod 2), (A38) 



7T 2 s 



TRIM re 



because we can evaluate 9 by writing the Hamiltonian as 
the k x — 7r cross-section of a Hamiltonian in one more 
dimension, and then determining the Chern number of 
that Hamiltonian^. The parity is given by (|A37I) . 



Appendix B: Tight Binding Models and Inversion 
Parity. 

A useful example to study is a tight-binding model 
with one atom per unit cell which has some orbitals that 
are even, |s,-), where i ranges from 1 to N even and some 
odd ones, \pi), where i continues from N even + l to N oc id + 



Let the Hamiltonian be 



H 



E 

i,i,R,A 



t v -(A) c l(R+A) Cj (R), (Bl) 



where Ci includes creation operators for both types of 
orbitals. 

A Bloch wave function ip a k (a is the band index) is 
determined by its values at the origin, <f>i = (si/pi\ip a k), 
for short) which are repeated in a plane-wave fashion to 
other sites. Inversion maps k > — k and <pi — > Ioij(f>j, 
where Ioij is a diagonal matrix of N even ones and N dd 
minus ones. 

There is a simple criterion for determining the parities 
of the occupied states at the TRIMs: if a state is odd, 
then all its s-wave components vanish. If it is even, all p- 
wave components vanish. In the generic case, it is enough 
to check a single s-wave component: 



T) a (K) = +1 if (si|V>are) ^ 
r} a (K) = -1 (probably) if (si|?/>are) = 0. 



(B2) 



If a component is exactly zero, it is reasonable to assume 
that it is because inversion symmetry forbids that com- 
ponent. 

The Hamiltonian used to illustrate the entanglement- 
rule in Fig. [5] is of this type. In momentum space Eq. 
([Bi"|) can be written £\ . cj(k)cj(k)/f i: ,-(k). The Hamil- 
tonian can be built from two subsystems Hi , H 2 , whose 
occupied states are given by the upper and lower signs in 
Fig. [2j Each Hamiltonian can be gapped because both 
sets of signs satisfy the constraint Y[ K ? ?re by itself. 

Each Hamiltonian H r can be constructed in a space 
with one even state \s r ) and one odd state \p r ). For the 
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Hamiltonian to be inversion symmetric it must commute 
with a z , so the diagonal matrix elements are even and 
the off-diagonal matrix elements are odd functions of k. 

The two Hamiltonians can be combined together by 
including a small mixing between the two systems. The 
final Hamiltonian has the form: 



H(k) 



H x (k) tl 
tt H 2 (k) 



(B3) 



where 



Hi = A k a z 
H 2 = D k a z 



B\,a x 



F k <Jy 



(B4) 



and 



Aic = (1 — cos (k x )) cos (k y ) + (1 + cos (k x )) cos (k y — k z ) 
£?k = (1 — cos (k x )) sin (k y ) + (1 + cos (k x )) sin (k y — k z ) 
\ = sin (k x ) ,E k = - sin (kx) , F k = sin (k z ) 



Q 



- cos (k x ) + cos (k z ) — 1. 



The parities of the occupied states at n of Hi and Hi un- 
der inversion in the point \x arc equal to the signs shown 
in Fig. [21 Choosing this point to define the inversion par- 
ities is necessary because the entanglement spectrum is 
obtained by cutting the system between two planes of the 
atoms. 



Appendix C: Monopoles 

Consider a Hamiltonian for which the a th and a + 1 st 
bands arc close to being degenerate at ko. Then all states 
except the two nearly degenerate states are unimportant, 
and the spectrum can be described by a 2 x 2 matrix 



H(k) 



E + A(k) B(k) 
B*(k) E + C(k) 



(CI) 



The eigenvalues are degenerate wherever A(k) = C(k) 
and B(k) = 0. This gives three equations in three vari- 
ables (since B is complex), so generically, a solution may 
be found in three dimensional space, as pointed out by 
Von Neumann and Wigner. 

Now at the location of the degeneracy (which we call 
kd) there is a magnetic monopole in the Berry magnetic 
fields B a and B a+ i of the two bands. The degeneracy 
point (or Dirac point) has a handedness 5 = ±1, which 
determines the charge of the monopoles. Bands a and 
a + 1 have opposite monopoles, of charge 



Qa = 2n6 

Qa+1 = -2nS. 



(C2) 



This monopole and antimonopole are glued together (at 
the point k^) but they do not annihilate because they 
are monopoles in different magnetic fields. 



To understand why the monopoles exist, write the 
Hamiltonian in terms of Pauli matrices, then expand the 
coefficients to linear order in k, H(k) as Eq + vo • k + 
52i=i( v i ' koi). The handedness is defined by 

S = sign[vi • (v 2 x v 3 )]. (C3) 

Define the oblique coordinate system Ki = (k — kd) ■ Vj. 
Now change to spherical coordinates, K% = K cos 9, Ki = 
K sin 9 cos <j>, K 2 — K sin 9 sin <f>. The eigenvalues are 
Eq + vo • k ± K. Note that if vo = 0, then the graph 
of the energy as a function of Ki , K 2 with K 3 set to 
is a cone. In accordance with the degrees-of-frccdom- 
argument given above, fixing K% at a random value (and 
thus reducing the number of parameters to 2) gives a 
hyperboloid, without any touching between the bands. 

The Berry magnetic field in band a is defined as B a = 
curlA a where A a is defined in Eq. (fTST) ; note that the 
distinction between u and ip does not matter for the pur- 
pose of calculating the magnetic field. The eigenvectors 

. / cos Ue-W \ , . /sin Ue~^ 

are ip a+ i = . 2 lo and ip a = i Q 

\ sin ^9 J \ — cos j9 

The Berry connections for these states are A a+1 (K) = 

\ cot |^ and A a (K) = ^ tan |0. These fields are familiar 

as the vector potentials for Dirac monopoles^- they have 

a flux-tube, or Dirac string, approaching the origin along 

the positive and negative K% axis. These monopoles have 

fluxes of (jij)B Q +i = — 2ir and (£j)B a = 2ir respectively. 

The calculation just completed is incorrect half the 
time- the magnetic flux can be calculated in any coor- 
dinate system provided it is right-handed; it is a pseu- 
doscalar. Hence, the result just obtained is correct when 
5 = +1. When 5 = — 1, the oblique coordinate system 
is left-handed, so the signs of the monopoles should be 
flipped, leading to Eq. (fC2|) . 

Now let us show that the curves defined by s a (k) = 
and s a+ i(k) = 0. (This appears in Sec. IIII B\ ) The 
function s a (k) is defined as (s|V> a k) = (s\ t)(t lV>ak) + 
(s| I) (| \tpak) (since only the two states | f) and | I) that 



are represented by 



in the effective theory 



are important near k^). So suppose 



(f \s) \ ( rns2 P -' ; /5 

a s) 



We then find that s a (k) = only if 9 = a and 4> = /3, 
and that s a+1 (k) = 0if^ = 7r-a and <j> = tt + /3. That is, 
the two curves are rays meeting at the Dirac point from 
opposite directions. 

Now let us see why a system cannot remain insulat- 
ing when n changes. When an even and odd state at a 
TRIM interchange, a pair of monopoles forms or annihi- 
lates. This may be seen using an effective Hamiltonian 
near the TRIM (say it is k = 0). Eq. (|C1[) can be recy- 
cled for this problem. Since IqH^Iq 1 = H(—k), and 

= ) the diagonal entries have to be even and the off- 
diagonal entries have to be odd functions of k. Assume 
Eq and the trace are equal to zero. (The trace just shifts 
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both bands by a smooth function and does not affect 
the degeneracy.) We wish to study how the dispersions 
change when A(0) — (7(0) = A, the difference between 
the two energies at the TRIMs, changes sign. 

Let A(k) = y + /(k) = — (7(k); / is a quadratic func- 
tion to lowest order. B(k) is a linear function of k, so 
we may choose a coordinate system where the real and 
imaginary parts are K\ and K 2 , B(k) = K\ — iK 2 . The 

dispersion is therefore ±v/0f + /(k)) 2 + K\ + K\. The 
states are degenerate if K = Kd where Kdi = Kd2 = 
and A = -2/(1^3,0,0) = -2aKj s , say. Thus, if 
a > 0, there are band-touchings when A is negative (with 

Kd3 = ± an d no crossings when A is positive. So 
if A changes from positive to negative, two monopolcs 
appear. The crossings move out from the origin to a 
distance proportional to the square root of A. 

Now reexpand the Hamiltonian to linear order around 
K<j. We find that it is equal to 01/^3(^3 — K^)cr z + 
K 1 a x + K 2 a y . (Notice that we have set Ki = K 2 = in 
f(Ki, K 2 , K 3 ) because we are interested in values of the 
if's such that K\,K 2 -C K3. The orders of magnitude 
such that all the terms under the square root in the dis- 
persion have the same magnitude are K\ ~ A, K 2 ~ A, 
K3 ~ a/A.) The two cone points have opposite hand- 
edness since has opposite signs. This was what we 
expected since the monopole charge has to be conserved. 

Appendix D: Entanglement Spectrum 

The flat band Hamiltonian is defined in terms of the 
correlation function C(xi,x 2 ) = (tp(x 2 Y '^(xi)) . (Once 
the results are derived in one dimension, the three di- 
mensional results can be written out by including factors 
of e ikj -' r .) 

Because the correlation function decays exponentially 
(for an insulator), it is reasonable to think of it as a 
hopping matrix for an electron system: 

[i?flat0](x) = ±<f>(x) - J C(x,x')<t>(x')dx'. (Dl) 

The eigenfunctions of this Hamiltonian are the same as 
the eigenfunctions of the physical system but the eigen- 
values are different. Each unoccupied state </> 7 has eigen- 
value i and each occupied state belongs to an eigenvalue 
-|, since (7(xi,x 2 ) = 2~] 7 ^ 7 (xi)0 7 (x 2 )*. 

When a surface at x = is introduced we can split 
the wave function into two parts, x > and x < 0, and 
represent them as the top and bottom halves of state- 
vectors. Then the correlation function has four parts, 

V <^RL ^R J 

The entanglement eigenstates /f turn out2£ to be the 
eigenfunctions of Hl = 5I — Cx- 



The eigenvalues are called \ — pi(k±). These eigen- 
values arc between ±| because pi is the probability that 
an electron occupies the state //'. Any eigenvalue in the 
range < p, < 1 is a state in the gap, so the wavefunction 
is localized near the surface. There are infinitely many 
surface states like this. In three dimensions, if they are 
graphed as functions of k^, the highest bands converge 
to±i. 

The entanglement energy spectrum is not quite the 
same as pt but is related by the transformation e = 
2tanh~ 1 (l — 2p) which sends the limiting energies to ±00; 
p = h corresponds to e = 0. The graphs of the entan- 
glement energies look more normal-they retain a finite 
spacing, although there are infinitely many of them. 

Since C has 1 and as eigenvalues, C 2 = C, giv- 
ing four matrix equations. For showing how to pair 
eigenstates^ of Cx and Cr, the relevant equation is 
Chl(1 — Cl) — CrCrl. Given an eigenfunction /£ of 
Cl with eigenvalue pn one can obtain an eigenvector of 
Cr with eigenvalue 1 — Pi via the transformation M 

/f (x) = [Mfhix) = . f ] = C RL {x,x')f*{x'). 

(D3) 

The prefactor is inserted to ensure that //* is normalized. 
Because of how p transforms, tLi = —£Ri ■ (One can check 
that M is a unitary transformation, which can be written 

in matrix form M = , » 1 » Crl ) ■ The eigenvalues of 

y Cr— C\ 

ft 4 and fP for H L and H R are ±{\ -pA and e Ll = -e m . 

Furthermore, F, = ^/p^/ ^ L + a/I — PifP is an occupied 
state because it satisfies C\Fi) = \Fj). (Eq. (|24|) restates 
this obscurely.) States of this form give a basis for all 
the occupied states. The ground state is therefore given 
by Yliiy/Pil} + x/ 1 ~Pi r \)\ vac ) wnc rc l],rj create the 
states /j L and fp respectively. Cross-multiplying gives 
the Schmidt decomposition, Eq. (|2"2"T) . The state Fj is 
definitely occupied by an electron. This electron is on 
the left with probability pi and on the right with prob- 
ability 1 — pi. A term in the Schmidt decomposition is 
obtained by making a choice, for each mode, whether the 
electron is on the left or right. 

If a state is filled on the left side, then the correspond- 
ing state must be left empty on the right. Thus the 
Schmidt states are \q)l = ELeA ^l\ vac )L an d \ch)r = 
Y[i£A r l\ vac ) R where A is a set of states. The weight of 
this state is H ieA Pi IL^ 1 ~ P*)- 

To reinterpret the fluctuations as statistical fluctua- 
tions of the system on the left, imagine covering the 
right half. Then electrons disappear when they cross 
the boundary. The factors of r\ correspond to holes in 
\q)l- The weight of a state can be written in terms 
of just the occupied states on the left by factoring out 
\ ■= 11,(1 - pi). The weight is then flieA T^i7 wn ich 
looks like a Boltzmann distribution if e~ tLi := , Pi , ex- 
plaining where the definition of en comes from. 

The maximum weight Schmidt state is obtained by 
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placing electrons on the left half when pi > \ and on 
the right half when pi < \. If we can only see the left 
half, this is equivalent to filling all the negative "energy" 
states. The negative energy states on the right are also 
filled since eju = — en- So the maximum weight state 
is the product \G)l\G)r of the ground states of Hl and 

Appendix E: Parity of Arcs through a TRIM 

Suppose AiV e (rej_) is prescribed at a certain TRIM. 
The number of arcs passing through this TRIM must be 
equal to Av modulo 2. To show this, we use the k ■ p 
effective Hamiltonian in the space of states that have en- 
ergy zero at k± to determine how the energies vary away 
from k±_ . Suppose for simplicity that all the modes at the 
TRIM have the same I e parity (the generic case). Then 
particle-hole symmetry implies that the effective Hamil- 
tonian is odd in kj_, to leading order H(k±) = A x k x + 
A y ky+. . . . Hence H(k±) = \k\(A x cos 9+A y sin 8), in po- 
lar coordinates. The energy-dependence has a cone-like 
structure: €i(k,9) = \k\fi(9) where fi(B) are the eigen- 
values of the periodic Hamiltonian Ac cos + A y sm9. 
Now the particle-hole symmetry implies that the bands 
come in pairs satisfying e,(k^) = — &£'(— kj_), or in other 
words that fi(9 + ir) = —fAN e +i-i(8)- Thus between 9 
and 0+7r, the energies must be turned upside down. This 
relates the dispersions fa in pairs, except for the middle 
one /^(Ajv c +i) (if AiV e is odd) which is related to itself. 
This mode changes sign from to it by the symmetry, so 
it crosses through at an odd number 2fc + l of values of 9 
in between, and crosses zero again at 2k + 1 points 180° 
away. When the solutions to ej(kj_) = are graphed, 
these crossings correspond to 2fc + l arcs passing through 
the TRIM. The other pairs of modes ft and f^+i-i to- 
gether give rise to an even number of arcs. Thus, the 
parity of the number of Fermi arcs is equal to the parity 
of AN e . 

Though the parity of the number of arcs is determined 
by AiV e 's parity, the precise number of arcs is not. For 
example, when AN e = 2, there may be no zero-energy 
states away from the TRIM, as in the Dirac equation. 

Appendix F: Polarization 

From the point of view of standard electrostatics, one 
expects a cubic crystal, aligned with the x, y, and z axes, 
with polarization P, to have a surface charge of P x per 
unit cell on the yz surfaces. However, in general, a ma- 
terial may have some stray charges on the surface. 

In some situations, the surface charge is expected to 
be determined (almost) by the theoretical value of the 
polarization™ This does not always happen because 
there may be extra charge trapped on the surface. But 
for a clean surface (that is, a perfectly periodic one) the 
ambiguity can be reduced: the charge density per unit 



cell is given by P x up to an integer multiple of e, if the 
surface is gapped. 

If there are no modes at the Fermi energy on the sur- 
face, then P x = Q x + ke (the surface charge per unit 
cell). The surface charge may be changed by ke by filling 
k surface bands. It is not possible to add a fractional mul- 
tiple of e to each cell, since then the surface will become 
conducting, according to the theorem that an insulator 
must have a whole number of electrons per unit cell, un- 
less strong interactions produce an unusual phase. (Since 
we are assuming time reversal to be broken, generically, 
an odd number of electrons can form an insulator.) 

This prediction for the surface charge may be general- 
ized to allow for a metallic surface. In this case there is a 
two dimensional Fermi surface describing the modes on 
the surface. The polarization is 

P* = Q x - ^% (mod e). (Fl) 

where Af s is the area of the surface arcs. The second 
term can be interpreted as the charge that needs to be 
removed to make the surface insulating. 

If a crystal has a fractional polarization of P x , then 
there are three possibilities. The surface may be electri- 
cally charged (with a density of P x + ke per unit cell for 
some integer fc), or it may be metallic. The surface may 
also reconstruct, so that there is a charge density wave. 
When P x is a simple fraction, this is very likely, since 
then the charge density wave would be commensurate. 
For example, when P x — | (as expected for inversion- 
symmetric insulators) the surface could have a period 
two reconstruction. The surface might also be metal- 
lic, but a big spontaneous surface charge seems unlikely. 
Since the "intrinsic" polarization may be exchanged for 
a surface property, a scanning tunneling microscope may 
be the best tool for observing it. 

The polarization may be determined by noting that the 
entanglement Hamiltonian has to satisfy the same 
constraints on its surface charge. We can just determine 
Q x and Af s for the ground state of H^. Assume Gr = 0. 
Then the number of arcs through each TRIM has the 
same parity, either even or odd (according to Eq. ([7])). 
The Fermi surface covers half the Brillouin zone if this 
number is odd, so according to Eq. ([231 . 

^ = i(AiV(0,0,0) + AA^,0,0)) 

71 ~ 

ee --P x (mod 1) (F2) 

where P x is defined by Eq. (fTS)) . The second line uses 
AN(k) =n- 2n (n) 

Now, if there are no nuclei at x = 0, then as in Sec. 
IIV1 one can use the neutrality of the Schmidt spectrum 
and symmetry to show that Q x = 0. If on x = there 
are nuclei of total atomic number per cell Z , imagine 
taking these nuclei out of the system. This leaves behind 
a charge of Z$e which must be divided evenly between the 
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left and right half. Focus on the left half of the system; 
it has a charge of Q x = 

Now substitute Q x and Af s into Eq. (|F1|) . Combine 
the n term in Eq. (jF2j) with to get (z °~" )e . By 
neutrality n is the total atomic number per unit cell, so 
n — Zq is the atomic number of the nuclei not on x = 0, 
which is congruent mod 2 to Zi (the number of nuclei 

at x = h) because the other nuclei come in pairs. Hence 
the polarization is eP e — ^Zie. This agrees with Eq. 
([5]) because the dipole moment of the nuclei relative to 
x = is — \Z\_e. Nuclei not on the two special planes 
and i come in pairs and cancel out. (We are defining the 
dipole moment of the nuclei to be the sum over a unit 
cell bounded by — ^ < x < 4 in the x-direction.) 

Appendix G: Frozen Crystals 

This appendix determines what subset of the space of 
n vectors is spanned by integer combinations of frozen 
crystals. Let f p be vectors corresponding to systems with 
a single fixed electron in each unit cell displace by p from 
the Bravais lattice. 

For each of the 7 nonzero polarizations p, f p is repre- 
sented by a vector with four ones and four zeros (f p (k) = 
^■p(mod 2)). This assumes the electron to be in an even 
orbital. For p = 0, take the electron to be in an odd or- 
bital instead so that fo(/s) = 1- 

The goal is now to determine what vectors are integer 
linear combinations of the f 's. There is a coordinate sys- 
tem for Z 8 where this is easy to solve One has to find a 
set of vectors Vi . . . Vs such that it is a basis for Z 8 , and 
also niVi, ri2V2, . . . , ngVg is a basis for the frozen vectors 
(where m, . . .ns are certain integers). Then if a vector 
n Q is represented by aivi + ■ ■ ■ + agVg in the new coor- 
dinate system, the criteria that it is a frozen vector are 
simple-a^ has to be a multiple of rij. The classification 
theorem for finitely generated abelian groups describes 
an algorithm for finding such bases. 

To find the basis, take an 8 x 8 matrix whose columns 
are the f's and do column operations. The only opera- 



tions that are allowed are adding or subtracting multiples 
of one column to another or changing the sign of a col- 
umn. These operations do not change the lattice spanned 
by the f's. (They can be inverted without dividing by in- 
tegers.) 

This process leads to the following basis: vo, v x , 2v xy , 
and 4,v xyz and vectors symmetric with these. Here,v 
is the vector with ones at all corners of the cube, v x is 
the vector with ones on the face of the cube defined by 
k-x = 7T (and zeros elsewhere), v xy is the vector with ones 
on the edge defined by n x = n y = tv, and ~v xyz is the 
vector with a one at the vertex 

The procedure for changing the basis from 
f ,fi,fij,%cyz to v ,Vi, 2Vij,4v xyz (where i,j run 
over x, y, z) consists of changing one of the four groups 
of f's to v's at a time. First fo = Vo,fj = v, so we 
can just rename the sequence of vectors Vo, Vj, £y, i xyz . 
Next, replace fi (s+s) by -(fi (ig+s) - v x - v y ), which is 
equal to 2~v xy . This is one of the column operations just 
described (it changes a vector in the third set by adding 
other vectors to it). Do the same for the other pairs of 
x,y,z. Last replace f|(a+g+5) by 

4v xyz = fl(x+ji+5) - Vx - Vy - V 2 + 2w xy + 2w yz + 2w xz . 

Each of these relations makes the ones cancel out except 
on just the right vertex or edge of the cube. 

Some thought leads to a similar demonstration that 
vo, Vjj and v xyz (without the factors of 2 and 4) span 
all combinations of 8 integers. Therefore any vector n Q 
can be decomposed as 

Oiij^ij + a xyz ~v xyz (Gl) 

i i<j 

where Ojj is an even integer and a xyz is a multiple of 
four. An even part of atj and a multiple of 4 contained 
in a xyz can be combined with the ao,ai-terms to form 
a vector that represents an inert insulator f , leaving the 
remainder given in Eq. (|21j) . 

The concise statement of this result is that the quotient 
of Z 8 by the span of the f's is x Z4. 
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This follows from the relation between topological insu- 
lators and the homotopy groups of Grassmann spaces 
Qn,m', the essential fact is that the homotopy group 
T2s (Gs,n-s) = 1 when N is sufficiently large—. 

27 There are technical problems when an orbital \a) is at a 
half-lattice vector-choosing which copy of it belongs to the 
unit cell breaks inversion symmetry. The simplest solution 
is to slice the orbital into two parts that are images of each 
other under inversion, \a) = -^=(|ai) ± |«2)). The opposite 
combination of these orbitals is then assigned a very large 
energy so that it has no effect. 

28 In a sequence such as U(1),U(2), etc., the homotopy 
groups of the first few elements are irregular, but even- 
tually stablizes farther in to the sequence. 

29 Suppose in particular that we are working with a hopping 
model consisting of a chain of identical sites each with mul- 
tiple orbitals. We have to make a cut on a bond midway 
between two sites, while defining inversion parity relative 
to sites is the most natural choice. Defining the parities 
this way, AN e = n (n) — n o (0) since the parities of states 
at k = 7r are switched by the change in inversion center 
from a bond to a site. This is equivalent to the result ob- 
tained in Sec. IIC of Refp— The factor of 2 in that article 
is present because the authors consider periodic boundary 
conditions, and cutting the chain creates two ends. 

3 In one dimension, this construction gives a basis of local- 
ized occupied states, somewhat like Wannier orbitals. Each 
of the f L and f R states, which have energies in the gap, 
are confined to some layer near the surface. Their sum is an 
eigenstate of //flat- (Unlike the Wannier basis, these states 
are all localized around the same bond.) The Hamiltonian 
//flat has localized eigenfunctions, unlike most Hamiltoni- 
ans, because the dispersion relation is flat, and so the group 
velocity is equal to zero. 

31 The Hamiltonian we are classifying (H'(k) = H(k)QH(0)) 
is one of the generalized Hamiltonians defined above; that 
is why a number of bands n e or n can be negative. 



